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FOREWORD 


This  work  was  performed  in  support  of  the  guided  projectile 
and  5"/54  ammunition  improvement  programs.  In  additon  to  the 
above  programs,  support  for  this  work  was  provided  by  the  Naval 
Ordnance  Systems  Command  under  ORDTASK  55A-501/090-1/UF  32-323-505. 

This  report  was  reviewed  by  Mr.,  D.  A.  Jones,  III,  Head  of  the 
Aeroballistics  Group,  Mr.  L.  M.  Williams,  III,  Head  of  the  Ballistics 
Division  and  Mr.  W.  R.  Chadwick,  Research  Aerodynamicist. 


-d 


1 


ABSTRACT 


Several  theoretical  and  empirical  methods  are  combined  into  a 
single  computer  program  to  predict  lift,  drag,  and  center  of 
pressure  on  bodies  of  revolution  at  subsonic,  transonic,  and 
supersonic  Mach  numbers.  The  body  geometries  can  be  quite  general 
in  that  pointed,  spherically  blunt,  or  truncated  noses  are  allowed 
as  well  as  discontinuities  in  nose  shape.  Particular  emphasis  is 
placed  on  methods  which  yield  accuracies  of  ninety  percent  or 
better  for  most  configurations  but  yet  are  comnutationally  fast. 
Theoretical  and  experimental  results  are  presented  for  several 
projectiles  and  a  computer  program  listing  is  included  as  an 
appendix. 
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INTRODUCTION 


In  the  past,  designers  have  relied  on  wind  tunnel  and  ballistic 
range  tests  to  predict  static  forces  and  moments  on  projectiles. 

This  is  not  only  very  expensive  but  also  quite  time  consuming  because 
of  the  man  hours  required  in  scheduling  and  performing  a  test  of  the 
above  nature.  At  most  test  facilities,  there  is  also  a  backlog  of 
work  of  about  three  to  six  months. 

It  is  believed  that  a  large  portion  of  wind  tunnel  and  ballistic 
range  tests  could  be  eliminated  (particularly  for  preliminary  and 
intermediate  design)  if  an  accurate  theoretical  method  were  available 
to  compute  static  forces  and  moments  throughout  the  Mach  number 
range.  More  important  though,  is  the  practical  use  of  such  a  method 
to  the  design  engineer  in  determining  the  configuration  which  is 
the  most  optimum  from  a  lift,  drag,  and  pitching  moment  standpoint  for 
his  given  design  goals.  Quite  often,  due  to  the  lack  of  such  a 
method  or  the  funds  for  wind  tunnel  testing,  less  than  optimum 
aerodynamic  configurations  are  used  to  accomplish  a  given  task. 

A  typical  example  is  the  external  shape  of  the  5"/38  projectile. 
According  to  the  work  of  reference  1,  the  range  of  that  projectile 
could  have  been  increased  by  more  than  fifty  percent  with  proper 
design. 

It  is  the  purpose  then  of  the  present  work  to  develop  a  general 
program  which  can  be  applied  to  the  body  of  the  guided  or  unguided 
projectile  to  predict  lift,  drag,  pitching  moment,  and  center  of 
pressure  over  the  Mach  number  range  of  current  interest,  0  <  M^  <  3. 
The  methods  used  in  the  development  of  the  program  should  be 
accurate  enough  to  replace  preliminary  and  intermediate  wind  tunnel 
tests  (accuracy  of  ninety  percent  or  better  for  most  configurations) 
but  yet  should  be  computationally  fast  enough  so  it  can  be  used  as 
an  efficient  design  tool. 

There  are  many  methods  available  in  any  particular  Mach  number 
region  to  compute  static  forces  and  moments  on  various  body  shapes. 
These  methods  range  in  complexity  from  exact  numerical  to  semi- 
empirical  and  the  body  shapes  vary  from  simple  pointed  cones  to 
complex  multi-stage  launch  vehicles.  However,  attempts  at  com¬ 
bining  the  various  methods  above  into  an  accurate  and  computationally 
fast  computer  program  have  been  scarce.  Saffell,  et  al2  developed 
a  method  for  predicting  static  aerodynamic  characteristics  for 
typical  missile  configurations  with  emphasis  placed  on  large  angles 
of  attack.  However,  the  drag  was  computed  by  handbook  techniques3 
and  slender  body  theory  was  used  for  the  lift  and  pitching  moment. 

As  a  result,  limited  accuracy  for  body  alone  aerodynamics  was 
obtained  using  this  method. 


Another  method  which  computes  forces  and  moments  throughout 
the  Mach  number  range  is  the  GE  "Spinner”  program4  designed 
specifically  for  projectiles.  This  program,  which  is  based  on 
empirical  correlations  as  a  function  of  nose  length,  boattail 
length,  and  overall  length,  gives  very  good  accuracy  for  most 
standard  shaped  projectiles.  However,  its  use  as  a  design  tool 
is  somewhat  limited  in  that  the  drag  of  a  given  length  nose  is 
the  same  no  matter  what  ogive  is  present  or  if  there  are  dis¬ 
continuities  present  along  the  nose.  The  same  statement  applies 
to  the  boattail  since  a  conical  boattail  of  from  5°  to  9°  is 
assumed  no  matter  what  the  boattail  shape  is.  Moreover,  no 
pressures  can  be  computed  by  the  GE  program  and  no  attempt  has 
been  made  to  include  nonlinear  angle  of  attack  effects. 

It  is  apparent  then,  from  the  above  discussion,  that  there  is 
a  definite  need  for  an  analytical  method  which  can  take  into 
account  nose  bluntness  and  ogive  shape,  disconti  nui  t j r?  i  long 

The  method  presented  herein  for  accomp'isningtnT^taSKre^es 
heavily  on  analytical  work  and  to  a  lesser  degree  on  empirical 
data,  As  such  it  is  believed  to  be  the  first  such  program  with 
major  emphasis  on  analytical  as  opposed  to  empirical  procedures. 

The  body  shapes  which  the  program  can  handle  should  be  general 
enough  so  that  most  projectile  and  missile  configurations  could 
be  handled  in  detail.  This  means  that  the  nose  may  be  pointed, 
truncated,  or  blunted  with  a  spherical  cap  and  that  the  nose  may 
have  two  ogives  present.  For  example,  on  a  typical  projectile 
the  fuze  has  one  contour  and  the  ogive  between  the  fuze  and 
shoulder  has  a  different  contour  with  a  discontinuity  in  between 
The  afterbody  should  consist  of  a  cylinder  followed  by  a  boattail 
or  flare.  A  typical  body  shape  along  with  the  coordinate  systems 
used  is  shown  in  Figure  1. 
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ANALYSIS 


A.  Wave  Drag 

Wave  drag  results  from  the  expansion  and  compression  of  the  air 
as  it  flows  over  the  body  surface.  Compression  of  the  air  is  seen 
in  the  form  of  shock  waves  which  first  occur  around  Mach  number  0.7 
to  0.9  depending  on  the  body  shape.  The  methods  used  to  calculate 
this  form  of  drag  differ  significantly  in  transonic  and  supersonic 
flow  and  thus  will  be  discussed  individually  below. 


Supersonic  Flow 

There  are  several  methods  available  for  calculating  the  suner- 
sonic  pressure  distribution  but  only  two  of  these  methods  hold 
promise  of  meeting  our  requirements  on  speed  of  computation  and 
accuracy  as  set  forth  in  the  introduction.  These  methods  are  the 
second  order  perturbation  theory  of  Van  Dyke5,6  and  the  second  order 
shock  expansion  theory7  modified  for  blunt  bodies  in  reference  8. 

Since  the  major  portion  of  the  flight  of  most  projectiles  is  in  the 
lower  supersonic  speed  regime  the  perturbation  approach  is  chosen 
because  it  is  more  accurate  than  shock  expansion  theory  at  these 
Mach  numbers.  However,  Van  Dyke's  theory  can  only  be  applied 
directly  to  bodies  where  the  slope  is  less  than  the  slope  of  the 
free-stream  Mach  lines.  Thus  for  blunt-nosed  configurations,  the 
perturbation  theory  is  combined  with  the  modified  Newtonian  Theory 
(the  means  for  combining  the  two  will  be  discussed  shortly). 

Before  discussing  the  combined  perturbation  Newtonian  approach 
a  brief  discussion  of  Van  Dyke's  theory  is  helpful. 

The  general  first  order  perturbation  problem  is:  (see  reference 
9  for  the  details  of  the  derivation) : 

TECHNICAL  LflfURT 
BLDG.  305 

4>rr  +  $r/r  +  4>09/r2  -  (M2  -1)  *xx  =  0  AJS3U»££tf  PROVISO  QRODSD,  un 

BIE1P-TL 


9  (0,r,9)  =  $x  (0 , r , 9)  =  0 


CD 


4>r  (x,R,0)  +  sin  a  cos  0  =  R'  [cos  a  +  4>x  (x,R,0)] 
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where  the  subscripts  indicate  partial  differentiation.  The  first 
order  problem  is  satisfied  exactly  by 


P  (x,r,e)  =  P  (x,r)  cos  a  +  ;  (x,r)  sin  a  cos  9  (2) 


where  the  first  term  corresponds  to  the  axial  flow  solution  and 
the  second  term  to  the  cross  flow  solution.  The  first  order  problem 
eq.  (1]  can  then  be  separated  into  an  axial  problem: 


<Prr  *  pT/r  -  &2PXX  =  0 


P  (0,r)  =  px  (o,r)  =  0 


(3) 


pr  (x,R)  =  R'  [1  -  p  (  X ,  R)  ] 


and  a  cross  flow  problem: 


4rr  +  ~  ^/r2  "  B^xx  =  0 

C  CO , r)  =  cx  (0,r)  =  0  (4) 

1  +  Cr  (x,R)  =  R'  rx  (X,R) 
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Without  going  into  the  details,  suffice  it  to  say  that  the  solutions 
of  eqs.  (3)  and  (4)  are  found  numerically  by  placing  a  distribution 
of  sources  and  doublets  respectively  along  the  x-axis. 

Van  Dyke  then  discovered  a  second  order  axial  solution  in  terms 
of  the  first  order  solution  ip.  Further,  since  disturbances  in  the 
cross  flow  plane  do  not  affect  the  pressure  as  much  as  disturbances 
in  the  axial  flow,  Van  Dyke  reasoned  that  it  would  be  quite  legiti¬ 
mate  physically  to  combine  this  second  order  axial  solution  with 
the  first  order  cross  flow  solution  of  Tsien1  •  to  form  a  hybrid  theory. 
Once  the  perturbation  velocities  ipx,  tpr  (second-order  axial)  ,  and  £x, 
and  £r  (first-order  crossflow)  are  computed  at  each  point  along  the 
body  surface  the  local  velocity  components  are: 


(cos  a)  (1  +  4ix)  +  (sin  a  cos  0)  (;x) 


(5a) 


v_  _  (cos  a)  (i|jr)  +  (sin  a  cos  0)  (1  +  cr) 
Vco 


(5b) 


w_  _  -  (sin  a  sin  8)  (1  +  ?/r) 


(5c) 


The  pressure  Coefficient  at  each  body  station  is  then: 


C  (x,0)  = 

P  ym£ 


1 


m2 


1 


U2  +  V2-t-W 


V2 


(6) 


Finally  the  force  coefficients  are: 


C 


A 


2 

TT  R£ 


Cp  (x,e)  de  dx 


(7) 


5 


c 


N 


1  |  Cp  (X,e)  cos  0  r  d0  dx 

*^0  *'0 


(8) 


C 


M 


,Z  rn 

Cp  (x,0)  cos  0  x  r  d0  dx 

/  I 

70  -o 


(9) 


and  the  center  of  pressure  in  calibers  from  the  nose  is 


"  -  CM  /  CN 


(10) 


It  should  be  pointed  out  that  in  the  actual  numerical  integration  of 
e1s-  (7),  (8)  and  (9)  the  integration  must  be  carried  out  in  segments 
of  the  body  between  each  discontinuity  due  to  the  discontinuous 
pressure  distribution. 

If  the  nose  is  pointed,  one  need  go  no  further.  But  if  the 
nose  is  truncated  or  is  blunted  with  a  spherical  cap  then  some  other 
method  must  be  used  to  determine  the  pressure  distribution  over  the 
truncated  portion.  The  method  used  herein  is  modified  Newtonian 
theory11.  Although  this  theory  is  derived  assuming  a  very  large 
Mach  number,  reasonable  values  for  the  pressure  coefficient  can  be 
obtained  over  a  portion  of  the  nose  even  at  low  supersonic  Mach 
numbers.  The  modified  Newtonian  pressure  coefficient  is 


sin' 


6 


(ID 


where  6  is  the  angle  between  a  tangent  to  the  local  body  surface 
and  the  freestream  direction  and  where  the  stagnation  pressure 
behind  a  normal  shock  is: 


6 


2  ( 

l"(Y+l)  Ml  1 

r  y+i  l 

Y  M2  ( 

L  2  J 

L  2yM|  -  (y-1)  -1 

(12) 


According  to  reference  12,  if  the  nose  is  truncated  then  the  pressure 
on  the  truncated  portion  is  only  about  ninety  percent  of  the  stagnation 
value  given  by  eq.  (12)  so  that  for  a  truncated  nose: 


0.9  C^1)  M“  ' 

2 


Y 

Y-l 


Y+l _ 

2y  Ml  -  (y-1) 


] 


(13) 


If  the  nose  has  a  spherical  cap  then  it  can  be  shown  that: 


6  =  sin  sin  B  cos  a  -  cos  8  cos  0  sin 


(14) 


where  tan  8  =  dr/dx. 

Then  combining  eqs.  (11)  and  (14)  one  obtains  for  a  spherical  nose 
cap: 


cp  (x»9)  =  Cpo  ^  sin2S  cos2a  -  sin2a  sing  cosB  cos0 
cos2S  cos 2 6  sin2a  ^ 


(15) 


where  Cp^  is  given  by  eq.  (12). 
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The  only  question  that  remains  now  so  far  as  the  supersonic 
Mach  number  region  is  concerned  is  where  does  the  modified  Newtonian 
theory  end  on  the  body  and  where  does  the  perturbation  theory  begin. 
To  determine  this  match  point  recall  that  the  slope  of  the  body 
surface  must  be  less  than  the  Mach  angle  to  apply  perturbation 
theory,  that  is 


6  <  sin 


(16) 


Thus,  the  upper  limit  of  the  perturbation  theory  is  5  =  sin  (l/'U) . 

Using  this  relation  in  eq.  (14)  and  assuming  a  spherical  nose  cap 
there  is  obtained  for  the  coordinates  of  the  point  below  which 
Newtonian  theory  must  be  applied: 


ru  = 


(V?£-i  cos  a  +  sin  a\ 


(17) 


x,,  = 


ru  tan 


a  +  rn  ( 


1  - 


1 


cos 


It  is  important  to  note  here  that  if  x  >  xu  Newtonian  theory  may 
still  be  applied  but  if  x  <  xu  perturbation  theory  cannot  be 
applied. 

The  limiting  angle  of  eq.  (16)  corresponding  to  the  coordinates 
of  eq.  (17)  is  shown  in  Figure  2  as  the  upper  curve.  Note  that  very 
large  cone  half  angles  can  be  computed  using  the  perturbation  theory 
at  the  lower  Mach  numbers.  However,  as  shown  by  Van  Dyke5  the  loss 
in  accuracy  of  perturbation  theory  increases  rapidly  as  the  angle  6 
is  increased.  Realistically,  since  at  an  angle  of  25°  -  30°  the 
error  is  still  slight  the  maximum  angle  6  for  which  perturbation 
theory  is  applied  should  not  exceed  these  values.  Based  on  these 
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accuracy  considerations,  the  Newtonian  theory  should  be  applied  for 
<5  values  outside  the  solid  line  boundary  of  Figure  2  and  perturba¬ 
tion  theory  within  the  boundary.  Now  the  match  point,  which  for 
the  present  work  will  be  defined  as  the  point  where  the  pressure 
coefficients  of  the  Newtonian  theory  and  the  perturbation  theory 
are  equal,  can  be  determined  as  the  solution  proceeds  downstream. 
For  body  stations  downstream  of  the  match  point,  perturbation 
pressures  are  used  in  the  force  coefficient  calculations  of  eqs . 
(7),  (8)  and  (9)  whereas  for  x  values  along  the  surface  less  than 
that  at  the  match  point  Newtonian  pressures  must  be  used. 


Transonic  Flow 

If  the  flow  is  transonic,  the  available  theories  for  the  wave 
drag  calculations  are  again  limited.  Here  the  main  limitations 
are  in  body  shape  because  there  does  not  appear  to  be  a  theoretical 
method  available  which  can  handle  the  blunted  nose  or  the  dis¬ 
continuities  along  the  body  surface.  Wu  and  Aoyoma13’11’  have 
developed  a  method  which  handles  tangent-ogive-cylinder-boattail 
configurations  at  zero  angle  of  attack  but  no  general  nose  geom¬ 
etries  can  be  used  as  is  the  case  in  supersonic  flow.  Thus  the 
approach  of  the  present  paper  will  be  to  calculate  the  wave  drag 
for  tangent  ogive  noses  of  various  lengths  throughout  the  transonic 
Mach  number  range  and  to  estimate  the  wave  drag  of  the  more  com¬ 
plicated  nose  geometry  based  on  these  results.  It  is  true  that 
the  accuracy  here  is  not  consistent  with  that  of  the  supersonic 
work  but  it  appears  from  the  results  (as  will  be  discussed  later) 
that  this  approach  is  justified,  at  least  for  noses  with  slight 
blunting  (rn/rb  <_  0.3)  . 

For  transonic  flow  the  perturbation  equation  (1)  has  an  addi¬ 
tional  term  so  that  for  a  =  0,  eq.  (1)  is  replaced  by: 


[l  -  M|  -  (y+ 1)  Mi,  #x  ]  4>xx  *  *r/ r  ♦  brr  =  0 


(18) 


Eq.  (18)  is  now  nonlinear  as  opposed  to  the  linear  eq.  (1)  used 
in  supersonic  flow.  Eq.  (18)  is  again  solved  numerically13  for 
the  velocity  potential  and  the  pressure  and  axial  force  coeffi¬ 
cients  calculated  by  eqs.  (6)  and  (7)  for  the  nose  of  the  nose¬ 
cy  linder-boattail  configuration.  Figure  3  gives  the  wave  drag 
obtained  by  solving  eq.  (18)  for  tangent  ogive  noses  of  various 
lengths  throughout  the  transonic  Mach  number  range.  For  a  given 
Mach  number  and  nose  length,  the  axial  force  coefficient  can  be 
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obtained  from  this  curve  by  interpolation.  If  the  pressure 
coefficients  along  the  body  surface  are  desired,  however,  the 
general  program  of  Wu  and  Aoyoma13  must  be  used. 

The  pressure  coefficient  on  the  boattail  at  zero  angle  of 
attack  in  transonic  flow1 11  is  given  by: 


,  1/2 

S  (X)  «  -  2  (X!  -  C)  r  _1  (X!  -C)?-  .  1  -  M?  -1  _  /dR>2 

5  /(Y+l)  M»2/3  L  25  (y+1)  Mcc2/3  (y+1)  M2  J  'dx' 


where  Xj  is  measured  from  the  shoulder  of  the  boattail  and 


2 

c2  =  25  (Y+l)  n£/3  ( I  A  _:__yL  +\  L  (  1  "  ) 

\  2  (Y+l)  M*  L4  V  (Y+l)  M2  / 


2/3 

4/3  _ 

1/2 

+  2  /  1  -  M 2  \ 

/  3  dR/dx  \ 

+  /  3  dR/dx  ' 

\ 

1 

M=o2/3  \(Y+1)  M^/ 

\  2  /y+ 1  / 

\  2Moo  A+l  , 

) 

) 

In  addition  to  the  restriction  of  zero  angle  of  attack,  eq.  (19) 
is  to  be  applied  for  1  <_  M  <1.2  (for  >_  1.2,  afterbody  wave 
drag  is  calculated  using  tffe  previous  supersonic  theory  for  the 
entire  body).  For  <  1,  experiment  shows  that  the  shock  first 
occurs  on  a  boattail  at  ~  0.95.  Accordingly,  wave  drag  will  be 
assumed  to  vary  linearly  from  zero  at  Moo  =  0.95  to  its  maximum 
value  at  ^  =  1.0  which  is  calculated  using  the  above  equation. 


(19) 
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B.  Skin  Friction  Drag 

The  boundary  layer  will  generally  be  turbulent  over  about  ninety 
percent  of  the  projectile  body  for  large  caliber  projectiles.  Since 
the  laminar  flow  region  is  usually  less  than  ten  percent  of  the 
total  surface  area,  it  will  be  assumed  the  entire  boundary  layer  is 
turbulent.  Under  this  assumption  the  total  or  mean  skin-friction 
coefficient,  Cfo0,  according  to  Van  Driest15  must  be  obtained  from: 


0-242  (sin 

A  (Cfji'2  (TW/TJ1/2 


C!  +  sin  -1  C2)  =  log1Q  (RNm  CfJ 


_^I_*_2ny°g10 


CTW/TJ 


(20) 


where  C.  =  2A  '  -  B _ 

(B2  +  4A2)1/2 


C2  =  B 

(B2  +  4A2)  ‘/2 


and 


A 


1/2 

(Y-l)  1 

2  VT„  ; 


B  =  1  *  (y-D/2  Nlj, 

Tw/Tro 


The  variable  n  of  eq.  (20)  is  the  power  in  the  power  viscosity  law: 


(21) 
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and  n  for  air  is  0.76.  Eq.  (20)  assumes  a  fully  developed  turbulent 
boundary  layer  with  zero  pressure  gradient  and  Prandtl  number  equal 
to  one. 

In  order  to  solve  eq.  (20)  for  the  mean  skin  friction  coefficient 
Cfa>,  one  must  have  values  for  Tw/Troj  RN  ,  and  M*,.  The  freestream 
Reynolds  number  is  simply 


RNoo  =  0t"  t 


To  relate  Tw/Too  to  the  freestream  Mach  number,  assume  the  wall  is 
adiabatic.  Defining  a  turbulent  recovery  factor  RT  by 


then 


(23) 


Too 


It  has  been  shown  that  the  recovery  factor  varies  as  the  cube  root 
of  the  Prandtl  number  (see  reference  16)  for  turbulent  flow  so  that: 


(24) 
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Recall  that  Van  Driest 's  Method  assumes  a  Prandtl  number  of  unity 
so  if  this  were  used  then  Rj  would  also  be  unity.  However,  the 
actual  value  of  Pr  =  0.73  so  that  the  previous  assumption  of 
Prandtl  number  one  can  be  compensated  for  somewhat  by  the  above 
recovery  factor  which  for  Pr  =  0.73  would  be  0.90.  Thus  eq.  (23) 
becomes : 


Tw/Too  -  1  +  0.9  III  ^ 
2 


(25) 


Then  for  a  given  set  of  freestream  conditions  (Ha,  ft*,,  Pco,  V*,)  one 
can  combine  eqs.  (22)  and  (25)  with  (20)  to  solve  for  Cf^.  The 
equation  must  be  solved  numerically  however,  since  Cf^  cannot  be 
solved  for  explicitly.  A  procedure  adaptable  to  equations  of  this 
type  is  the  well  known  Newton-Raphson  method  discussed  in  reference 
17. 


Once  the  mean  skin  friction  coefficient  has  been  determined  for 
a  given  set  of  freestream  conditions,  the  viscous  axial  force 
coefficient  is  simply: 


(26) 


The  wetted  area  Sw  is  the  total  surface  area  of  the  body  which  can 
be  integrated  numerically  given  a  set  of  body  coordinates. 


C.  Base  Drag 

Much  theoretical  work  has  been  performed  to  predict  base  pressure 
(references  18  -  22).  There  is  still  no  satisfactory  theory  available, 
however,  and  the  standard  practice  has  been  to  use  empirical  methods. 
This  is  the  approach  taken  here.  Figure  4  is  a  mean  curve  of  experi¬ 
mental  data  from  references  18,  19  and  23  -  29.  This  data  assumes  a 
long  cylindrical  afterbody  with  fully  developed  turbulent  boundary 
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layer  ahead  of  the  base.  There  could  be  deviations  from  this  curve 
due  to  low  body  fineness  ratio,  boattails,  angle  of  attack,  Reynolds 
number  and  surface  temperature.  Each  of  these  effects  will  be  dis¬ 
cussed  below. 

The  minimum  length  of  most  projectiles  is  about  four  calibers. 
According  to  references  23  and  29  the  base  pressure  at  low  super¬ 
sonic  Mach  numbers  is  essentially  unaffected  by  changes  in  body 
length  if  the  fineness  ratio  is  greater  than  four.  This  is  not 
true  at  high  supersonic  and  hypersonic  Mach  numbers  as  shown  by 
Love18.  But  since  the  main  interest  is  for  Mo,  3  the  effect  of 
overall  fineness  ratio  on  base  pressure  can  be  neglected. 

In  addition  to  the  above,  Love  shows  that  the  nose  shape  has 
little  effect  on  base  pressure  for  high  fineness  ratio  bodies. 

Thus,  for  bodies  of  fineness  ratio  of  four  or  greater  the  effect 
of  nose  shape  and  total  length  on  base  pressure  can  be  neglected. 

The  base  pressure  is  significantly  altered  by  the  presence 
of  a  boattail  so  that  this  change  must  be  accounted  for.  Probably 
the  most  simple  method  to  do  this  is  an  empirical  equation  given 
by  Stoney20, 


3 

caBa  =  'CPBA  (^) 


Eq.  (27)  can  be  used  throughout  the  entire  Mach  number  range  where 
Cp  is  the  base  pressure  given  by  the  curve  of  Figure  4.  An  alter¬ 
native  to  this  procedure  is  to  find  the  base  pressure  as  a 
function  of  boattail  angle  and  then  the  diameter  of  the  base  would 
be  squared  instead  of  cubed  as  in  equation  (27).  That  is 


Caba 
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where  C£  is  the  base  pressure  coefficient  for  a  given  boattail 
angle,  ihis  requires  knowing  CpB^  however  which  is  not  always 
available.  Because  of  this,  eq.  (27)  will  be  used. 

It  has  been  shown  in  many  works21’30  that  the  base  pressure 
is  essentially  independent  of  Reynolds'  numbers,  R^j,  if  the  boundary 
layer  ahead  of  the  base  is  fully  developed  turbulent  flow.  A 
turbulent  boundary  layer  usually  occurs  for  R\j  of  500,000  to  750,000 
depending  on  the  roughness  of  the  body  surface.  The  minimum  R^ 
ahead  of  the  base  one  would  expect  to  encounter  on  the  present 
bodies  would  be  about  1,000,000.  Moreover,  most  projectiles  have 
various  intrusions  and  protrusions  such  as  on  a  fuze  which  tends 
to  promote  boundary  layer  separation.  In  view  of  these  practical 
considerations,  Reynolds  number  effects  on  base  pressure  may  safely 
be  neglected. 

The  same  arguments  as  the  ones  above  hold  for  surface  tempera¬ 
ture  as  well.  Thus  in  addition  to  Reynolds  number  effects,  surface 
temperature  effects  on  base  pressure  need  not  be  accounted  for. 

The  effect  of  angle  of  attack  on  base  pressure  is  to  lower  the 
base  pressure  and  hence  to  increase  the  base  drag.  For  bodies 
without  fins,  the  amount  of  this  decrease  is  dependent  mainly  on 
freestream  Mach  number.  If  a  is  given  in  degrees  then  an  empirical 
relation  for  the  change  in  base  pressure  coefficient  due  to  angle 
of  attack  is  given  by 


[ACpJ  =  -(.012  -  .0036M*,)  a 

a 


(29) 


Eq.  (29)  was  derived  from  a  compilation  of  experimental  data  pre¬ 
sented  in  Figures  7  through  15  of  reference  23.  The  base  drag 
coefficient  thus  becomes,  in  light  of  eqs.  (27)  and  (29): 
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(.012  -  . 0036Muo)a 


(30) 
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D.  Viscous  Separation  and  Rotating  Band  Drag 

Figure  5A  is  a  plot  of  forebody  drag  coefficient  as  a  function 
of  cone  half  angle  from  data  taken  from  reference  31.  Since  the 
skin  friction  drag  coefficient  is  about  0.02  for  this  case,  it  can 
be  subtracted  from  the  curve  of  Figure  5A  to  yield  the  pressure 
drag  coefficient.  Note  that  the  freestream  Mach  number  is  0.4, 
low  enough  so  that  no  appreciable  compressibility  effects  occur. 

The  question  therefore  arises  as  to  the  origin  of  this  type  of 
drag,  since  it  is  not  compressibility  or  skin  friction  drag.  It 
is  in  fact  viscous  separation  drag.  For  very  large  cone  half  angles, 
0C,  the  flow  over  the  cone,  instead  of  remaining  attached,  separates 
due  to  the  very  strong  adverse  pressure  gradient  and  reattaches 
downstream.  This  separation  prevents  the  pressure  from  decreasing 
as  much  as  it  would  in  inviscid  flow  and  produces  a  drag.  Oddly 
enough,  this  phenomenon  does  not  occur  on  ogives  or  on  spherical 
surfaces,  apparently  due  to  body  curvature  effects  on  the  boundary 
layer.  As  a  result,  one  can  derive  an  empirical  expression  for 
this  viscous  separation  drag,  where  the  important  parameter  is  the 
angle  6*  which  the  nose  makes  with  the  shoulder  of  the  afterbody. 
Based  on  Figure  5A  this  relation  is 


CAv.s  =  .012  (6*  -  10°)  ;  6*  >  10° 


=  0;  6*  <  10° 


(31) 


with  5*  in  degrees  and  with  6*  =  9C  for  a  conical  nose. 

Reference  1  gives  the  measured  effect  of  a  rotating  band  on  drag. 

The  particular  rotating  band  used  in  those  wind  tunnel  tests  had  a 
mean  height  of  about  0.024  calibers.  An  expression  which  functionaries 
the  above  results  for  drag  increment  due  to  a  rotating  band  is  given 
by: 


cArb  -  (ACa)  (H)/.01 


(32) 
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where  H  is  the  mean  height  of  the  band  in  calibers  and  AC^  is  the 
increment  in  axial  force  for  an  H  of  0.01  caliber  given  in  Figure 
5B.  Although  the  eq.  (32)  was  derived  for  a  particular  band,  it 
checks  well  with  the  results  of  Charters32  for  a  different  band 
geometry. 


E.  Inviscid  Lifting  Properties 

At  supersonic  Mach  numbers  the  inviscid  lift,  pitching  moment, 
and  center  of  pressure  are  calculated  using  Tsien's  first  order 
cross  flow  theory  which  was  discussed  earlier  in  conjunction  with 
Van  Dyke's  second  order  axial  solution.  This  method  is  adequate 
for  small  angles  of  attack  where  viscous  effects  are  negligible. 

At  subsonic  and  transonic  Mach  numbers  the  lifting  properties 
are  more  difficult  to  obtain.  For  subsonic  velocities  the  lift 
could  be  calculated  by  perturbation  theory33  but  since  projectiles 
rarely  fly  at  Mach  numbers  less  than  0.7,  a  formulation  on  this 
basis  was  not  justified.  An  alternative  would  be  slender  body 
theory  but  the  accuracy  of  this  approach  is  inadequate.  In  light 
of  the  above  reasoning,  a  semi-empirical  method  for  normal  force 
characteristics  was  derived  based  on  nose  length,  afterbody  length, 
and  boattail  shape.  This  method  was  then  extended  through  the 
transonic  Mach  number  range  since  the  state-of-the-art  in  transonic 
flow  does  not  allow  one  to  handle  the  general  body  shapes  or  flow 
conditions „ 

The  total  inviscid  normal  force  acting  on  the  body  may  be 
written 


Cc.v> 


(33) 


where  the  subscripts  n,  a,  and  B  stand  for  nose,  afterbody,  and 
boattail  respectively.  The  first  term  of  eq.  (33)  can  be  approxi¬ 
mated  by 


=  Cj  tan  5*  ♦  C2 


(34) 
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where  and  C2  are  given  in  Figure  6  as  a  function  of  Mach  number. 
This  relationship  was  determined  empirically  from  the  cone  results 
of  0wen631.  It  is  approximately  correct  for  SLn  >  1.5,  cone  blunt¬ 
ness  up  to  0.5,  and  Moo  <  1.2.  Note  that  the  angle  6*  in  eq.  (34) 
is  the  same  as  that  discussed  previously  in  eq.  (31). 

The  normal  force  coefficients  of  the  afterbody  and  boattail  can 
be  obtained  from  Figures  7  and  8  respectively.  Figure  7  was  derived 
analytically  in  the  transonic  Mach  range  from  the  method  of  Wu  and 
Aoyoma13  and  in  subsonic  flow  from  the  experimental  data  of  Spring31* 
and  Gwin33.  In  the  work  of  Spring  and  Gwin  above,  the  normal  force 
of  the  nose  plus  afterbody  was  given  but  the  nose  component  can  be 
subtracted  off  by  the  use  of  eq.  (34) .  The  boattail  normal  force 
coefficient  was  given  by  Washington36  but  he  stated  that  there  was 
not  enough  data  available  in  subsonic  and  transonic  flow.  Hence 
the  data  of  Washington  was  supplemented  by  the  175mm  Army  projectile37 
and  Improved  5"/54  Navy  projectile38  data  to  derive  the  general  curve 
of  Figure  8. 

Although  slender  body  theory  may  not  be  adequate  for  predicting 
the  normal  force  coefficient  it  appears  to  predict  the  center  of 
pressure  of  the  nose  and  boattail  lift  components  quite  adequately. 
According  to  slender  body  theory  the  center  of  pressure  of  the 
nose  is 
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or 


(xcp)B  - 


£  .  (Vol)B 


IT  R‘ 


(36) 


The  center  of  pressure  of  the  afterbody  normal  force  was  calculated 
analytically  by  the  method  of  Wu  and  Aoyoma  in  transonic  flow  and 
assumed  to  have  the  same  value  in  subsonic  flow.  Figure  9  is  a  plot 
°f  (xcp^a^a  versus  afterbody  length  measured  at  the  point  where 
the  afterbody  begins.  Mow  knowing  the  individual  lift  components 
and  their  center  of  pressure  locations,  one  can  compute  the  pitching 
moment  about  the  nose  as: 


(cN  ) 
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F.  Viscous  Lifting  Properties 

Strictly  speaking,  the  previous  discussion  on  inviscid  lifting 
properties  gave  CN  and  at  a  =  0  only.  If  a  >  0  then  there  is 
a  nonlinear  contribution  to  lift  and  hence  pitching  moment  due  to 
the  viscous  crossflow  of  velocity  V  =  Vro  sin  a.  Allen  and  Perkins39 
list  these  contributions  as: 


(ACN)vis  =  n  cdc  a2 

Sr 


(38) 
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(39) 


fACM)v.s  .  -r,  tdc(!E)(xp)  o' 


where  n  and  cd  are  given  in  Figure  10.  Note  that  the  cross  flow 
drag  coefficient  is  here  taken  to  be  a  function  of  Mach  number  only 
and  the  cross  flow  Reynolds  number  dependence  is  not  accounted  for. 
The  center  of  pressure  of  the  entire  configuration  should  then  be: 


xcp 


Cm  +  (ACm) vis 

CN  +  (ACN)vis 


(40) 


G.  Summary 

Figure  11  gives  a  summary  of  the  various  methods  used  in  each 
particular  Mach  number  region  to  compute  the  static  aerodynamics. 
As  may  be  seen,  major  emphasis  has  been  placed  on  analytical  as 
opposed  to  empirical  procedures. 
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RESULTS  AND  DISCUSSION 


A.  Numerical  Solutions 

A  computer  program  was  written  in  Fortran  IV  for  the  CDC  6700 
computer  to  solve  the  various  equations  discussed  in  the  analysis 
section  by  numerical  means.  The  various  methods  used  for  each 
individual  equation  are  the  same  as  those  discussed  in  the  refer¬ 
ences  pertaining  to  the  particular  equation  and  will  not  be 
repeated  here.  However,  mention  should  be  made  of  the  fact  that 
the  step  size  used  in  the  hybrid  theory  of  Van  Dyke  was  considerably 
smaller  than  he  suggested,  particularly  for  a  blunt  nosed  body  or 
behind  a  discontinuity.  For  example,  for  the  most  complicated  body 
shapes  as  many  as  200  points  were  placed  along  the  body  surface. 

Also  slight  oscillations  in  the  second  order  solution  were  found 
behind  a  corner  although  Van  Dyke  does  not  mention  these  details. 

Quite  often,  it  was  necessary  to  evaluate  an  integral  numerically 
or  to  compute  the  value  of  a  function  and  its  derivative  at  a  given 
point.  The  integration  was  carried  out  using  Simpson's  rule;  the 
interpolation  and  differentiation  using  a  five  point  Lagrange  scheme*7. 
Both  methods  have  truncation  errors  which  are  consistent  with  the 
accuracy  of  the  governing  set  of  flow  field  equations . 

The  computational  times  depend  on  how  complicated  the  body 
shapes  are  and  the  particular  Mach  number  of  interest.  The  longest 
computational  time  for  the  most  general  body  shape  computed  was 
less  than  half  a  minute  for  one  Mach  number.  For  most  configurations 
the  average  time  is  about  fifteen  seconds  per  Mach  number  for  Mo.  > 

1.2  and  about  five  seconds  per  Mach  number  for  Moo  <  1.2.  This 
assumes  of  course  that  a  table  look-up  procedure  is  used  in  the 
transonic  region  where  the  curves  of  Figure  3  are  input  as  data 
sets  as  opposed  to  solving  the  nonlinear  partial  differential 
equation  (18)  for  each  Mach  number.  If  the  aerodynamic  coefficients 
of  a  given  configuration  are  desired  throughout  the  entire  Mach 
number  range,  an  average  execution  time  of  two  minutes  is  required 
for  most  configurations  (ten  Mach  numbers) . 

A  detailed  discussion  of  the  computer  program  is  included  as 
Appendix  A.  The  various  input  and  output  parameters  are  defined 
and  a  listing  of  the  program  along  with  a  sample  output  are  also 
included  for  the  reader's  convenience. 


B.  Comparison  with  Experiment 

The  only  new  method  presented  in  the  current  work  is  the  com¬ 
bined  perturbation  -  Newtonian  theory  for  blunt  bodies,  It  is 
thus  of  interest  to  see  how  the  pressure  coefficients  along  the 
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surface  compare  with  experimental  data,  Figures  12  and  13  present 
two  typical  comparisons  at  H*,  =  1.5  and  2.96.  The  experimental 
data  is  taken  from  reference  8  which  combined  modified  Newtonian 
theory  with  shock  expansion  theory  to  compute  forces  on  blunted 
cones.  The  asymptotes  of  the  pressure  coefficient  in  each  of  the 
planes  computed  by  the  method  of  reference  8  are  also  indicated 
on  the  figures.  As  seen  in  the  figures  the  present  theory  predicts 
the  aerodynamics  much  better  than  shock  expansion  theory  at  Moo  = 

1.5  and  is  about  the  same  as  the  shock  expansion  approach  at  Moo  = 
2.96.  The  reason  for  this  is  that  the  basic  perturbation  theory 
was  derived  assuming  shock  free  flow  with  entropy  changes  slight; 
hence  the  theory  should  be  most  accurate  in  the  lower  supersonic 
speed  regime.  On  the  other  hand,  shock  expansion  theory  was 
derived  assuming  a  shock  present  and  so  one  would  expect  this 
method  to  be  better  than  perturbation  theory  as  M^  is  increased. 
Apparently,  the  crossover  point  is  around  Moo  =  2.5  to  3.0  so  that 
for  the  major  portion  of  the  supersonic  speed  range  of  interest 
in  the  present  analysis,  perturbation  theory  is  more  accurate. 

Another  interesting  point  in  Figure  12  is  the  discontinuity 
in  slope  of  the  pressure  coefficient  curve  which  occurs  at  the 
match  point.  This  is  because  in  the  expansion  region  on  the 
spherical  nose  the  perturbation  pressure  decreases  much  more  rapidly 
than  the  Newtonian  theory  and  as  a  result  the  overexpansion  region, 
which  occurs  at  low  supersonic  Mach  numbers,  is  accounted  for  quite 
well.  Note  that  the  match  point  is  different  in  each  plane  around 
the  surface  (x  =  0.11  to  0.14). 

One  of  the  questions  which  arises  in  the  development  of  a  general 
prediction  method  pertains  to  accuracy.  To  answer  this  question, 
force  coefficients  for  several  cases  were  computed  embracing 
variations  in  nose  bluntness,  Mach  number,  angle  of  attack, 
nose  length,  and  afterbody  length.  These  cases  are  presented  in 
Figures  14  through  19  along  with  experimental  data. 

The  first  of  these  cases  (Figure  14)  gives  the  axial  force 
coefficient,  normal  force  coefficient  derivative,  and  pitching 
moment  coefficient  derivative  as  a  function  of  nose  bluntness  for 
a  simple  blunted  cone  configuration.  Note  that  the  axial  force 
coefficient  includes  only  the  wave  plus  skin  friction  components 
because  the  base  drag  was  subtracted  out  of  the  given  set  of 
experimental  data.  An  important  point  here  is  that  very  good 
accuracy  is  obtained-  even  for  large  bluntness  ratios.  For  example, 
with  bluntness  rn/rg  =  0.6,  the  force  coefficients  are  in  error  by 
less  than  fifteen  percent  (this  is  assuming  of  course  there  is  no 
error  associated  with  the  experimental  data  which  is  not  exactly 
correct).  This  tends  to  verify  that  a  combined  perturbation- 
Newtonian  theory  can  be  used  successfully  for  blunt  configurations 
even  at  low  supersonic  Mach  numbers. 


22 


The  next  two  figures.  Figures  15  and  16,  compare  the  theoretical 
static  aerodynamic  coefficients  with  experiment  as  a  function  of 
Mach  number  for  blunted  cones  with  bluntness  ratios  of  0.2  and  0,406 
respectively.  Also  included  in  Figure  15  is  the  slender  body  theory. 
As  seen  by  the  error  comparisons  at  the  lower  part  of  Figure  15, 
accuracies  of  better  than  90  percent  can  be  obtained  throughout 
the  supersonic  Mach  number  range  for  the  force  coefficients.  Figure 
16  gives  the  aerodynamic  data  throughout  the  Mach  number  range  of 
interest.  Again  the  comparison  is  favorable  even  though  the  tran¬ 
sonic  wave  drag  was  computed  for  a  tangent  ogive  having  a  length 
equal  to  that  of  the  blunted  cone.  The  bluntness  causes  the  tran¬ 
sonic  drag  rise  to  start  at  a  lower  Mach  number  and  to  be  less 
abrupt  than  for  the  pointed  tangent  ogive. 

The  third  vairable  of  interest  is  angle  of  attack.  Figure  17 
presents  the  results  for  a  tangent  ogive  cylinder  of  nose  length 
six  calibers  and  total  length  fourteen  calibers.  Two  Mach  numbers 
are  considered,  =  1.5  and  =  2.5,  Again  the  results  are  quite 
good,  except  at  very  large  angles  of  attack. 

Figure  18  compares  the  force  coefficients  of  the  present  theory 
with  experiment  for  a  pointed  cone  of  various  lengths.  Also  shown 
for  comparison  with  the  M^  =  1.5  case  is  the  slender  body  theory. 
Although  perturbation  theory  is  usually  associated  with  nose  slender¬ 
ness  ratios  of  two  and  greater,  it  may,  nevertheless,  be  seen  that 
fair  accuracy  is  obtained  for  lengths  as  low  as  one.  This  corresponds 
to  a  cone  half  angle  of  about  twenty-five  degrees  which  is  the 
limiting  angle  used  in  the  combined  perturbation  -  Newtonian  theory 
as  shown  in  Figure  2.  For  the  =  0,5  case  eq,  (31 j  is  used  to 
calculate  the  viscous  separation  drag  which  is  added  to  the  skin 
friction  drag  to  get  the  total  forebody  drag  coefficient-  Using  this 
simple  formula,  excellent  agreement  with  experimental  data  is 
obtained 

The  final  variable  of  interest,  afterbody  length,  is  examined 
in  Figure  19.  The  nose  of  the  body  is  a  2.83  caliber  tangent 
ogive  For  zero  afterbody  length,  the  theory  agrees  with  experiment 
very  well.  However,  as  the  afterbody  length  increases  the  theory 
underestimates  the  afterbody  lift  at  the  lower  supersonic  Mach 
numbers  for  short  afterbody  lengths  and  at  the  higher  Mach  numbers 
for  long  afterbody  lengths.  This  loss  in  lift  predicted  by  the 
inviscid  theory  was  also  found  by  Buford 1,0  and  he  attributed  it  to 
boundary  layer  displacement  effects.  Even  so,  the  present  theory 
is  superior  to  slender  body  theory  which  gives  zero  lift  due  to  an 
afterbody. 

To  summarize  the  previous  five  figures,  one  could  say  in  general 
that  accuracies  of  ninety  percent  or  better  can  be  obtained  for  force 
coefficients  of  most  configurations.  However,  for  extreme  cases. 
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such  as  very  large  nose  bluntness  or  angle  of  attack,  the  accuracy 
will  be  decreased  and  the  amount  of  this  decrease  can  be  approxi¬ 
mated  from  Figures  14  through  19, 

The  next  several  figures  compare  theory  with  experiment  for 
several  spin  stabilized  projectiles.  Figures  20,  21  and  22  are 
Navy  projectiles:  the  5"/38  RAP  (Rocket  Assisted  Projectile)4', 
the  5"/54  projectile  42,  and  the  improved  5"/54  projecti le? B , 
which  has  a  longer  nose  and  boattail  than  the  standard  5"/ 54. 

Figures  23  and  24  are  Army  shapes:  the  175mm37  and  155mm43  projec¬ 
tiles  respectively.  For  the  detailed  drawings  and  other  aerodynamics 
of  these  shapes  the  reader  is  referred  to  the  references  cited  above. 

The  theoretical  zero  lift  drag  curve  of  the  5"/38  RaP  projectile 
along  with  three  sets  of  experimental  data’  and  an  Nlv'L  empirically 
derived  curve  are  shown  in  Figure  20.  Note  that  the  experimental 
data  varies  by  about  thirty  percent  for  Ho  <  1  and  by  ten  percent 
for  Ho  >  1.  The  theoretical  curve  tends  to  support  the  BRL  data 
subsonically  and  the  NOL  and  NWC  data  supersonically.  The  numbers 
in  parenthesis  are  the  factors  by  which  the  drag  curves  must  be 
multiplied  throughout  the  flight  of  the  projectile  to  match  actual 
range  firings.  The  NWL  empirical  curve  is  the  curve  which  is 
actually  used  in  range  predictions  due  to  the  failure  of  experimental 
data  to  predict  an  adequate  drag  curve.  This  empirical  curve  was 
derived  from  actual  range  firings.  It  should  be, therefore,  slightly 
high  because  of  yaw  induced  effects.  The  important  point  here  is 
that  for  this  particular  shell,  the  theory  agrees  better  with 
actual  range  firings  than  any  of  the  sets  of  experimental  data. 

Figures  21  and  22  give  the  static  aerodynamic  coefficients  for 
the  5"/S4  RAP  and  the  improved  5"/54  projectiles.  The  5"/54  RAP 
has  a  nose  length  of  about  2.5  calibers  and  a  boattail  of  0-5 
calibers  whereas  the  improved  round  has  a  2.75  caliber  nose  and  a 
1.0  caliber  boattail.  Also  the  5"/ 54  RAP  has  a  rotating  band 
whereas  the  other  shell  does  not.  For  both  shells,  excellent 
agreement  with  experimental  data  is  obtained  for  the  drag  coeffi¬ 
cient  throughout  the  entire  Mach  number  range.  Fair  agreement  is 
obtained  for  normal  force  coefficient  and  hence  pitching  moment 
and  center  of  pressure.  The  comparison  for  the  lifting  properties 
is  Mach  number  dependent:  in  the  low  supersonic  region  the  theory 
is  consistently  about  ten  percent  low  on  normal  force  whereas  at 
high  supersonic  speeds  it  compares  very  well  with  experiment. 

The  reason,  as  already  mentioned,  is  the  failure  of  the  inviscid 
theory  to  predict  afterbody  lift  correctly  at  low  supersonic  Mach 
numbers.  At  subsonic  and  transonic  Mach  numbers,  the  theory  does 
about  as  well  as  could  be  expected  considering  that  there  was  a 
considerable  amount  of  empirical  work  in  that  region. 
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For  boattailed  configurations,  such  as  the  5"/54  RAP  and  the 
Improved  5"/54,  it  was  found  necessary  to  account  approximately 
for  the  thick  boundary  layer  on  the  boattail.  This  was  done  by 
viewing  the  unpublished  shadow  graphs  obtained  in  conjunction 
with  the  work  of  reference  38.  Apparently,  a  maximum  boattail 
angle  of  six  degrees  can  be  allowed  before  boundary  layer  separa¬ 
tion  takes  place.  In  addition,  the  boundary  layer  displacement 
thickness  accounts  for  another  about  1/4  -  1/2  degree  decrease  in 
the  effective  boattail  angle.  These  two  results  were  used  to 
determine  effective  boattail  angles  on  all  boattailed  configura¬ 
tions.  Without  this  approximate  accounting  of  the  boundary  layer 
effect  on  the  boattail  shape,  the  lifting  properties  would  have 
been  in  error  by  an  additional  ten  percent  for  boattailed  con¬ 
figurations  . 

The  final  two  shells,  the  175  and  155mm,  are  considered  in 
Figures  23  and  24.  Again,  excellent  drag  predictions  are  made  by 
the  theory  and  good  predictions  are  made  for  normal  force  and 
center  of  pressure.  Intuitively,  one  would  expect  the  axial 
force  to  agree  better  with  experiment  than  the  lift  because  a 
second  order  approach  is  used  m  sunersonic  flow  for  tne  axial 
forces  whereas  a  first  order  cross-flow  theory  is  used  for  the 
normal  forces. 

Figure  25  presents  theoretical  results  for  the  five- inch 
guided  projectile.  The  nose  is  about  sixty  percent  blunt  with 
two  different  ogive  sections.  The  overall  length  is  10.58 
calibers  with  a  0.66  caliber  boattail,  7.24  caliber  afterbody 
and  2.68  caliber  nose.  Although  no  experimental  data  is  currently 
available  for  this  extreme  case,  it  is  expected  that  the  theory 
is  accurate  to  within  ten  percent  on  axial  force  and  twenty 
percent  on  lifting  properties. 
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CONCLUSIONS 


1.  A  general  method  has  been  developed  consisting  of  several  theore¬ 
tical  and  empirical  procedures  to  calculate  lift,  drag,  and  pitching 
moment  on  bodies  of  revolution  from  Mach  number  zero  to  about  three 
and  for  angles  of  attack  to  about  twenty  degrees. 

2.  Comparison  of  this  method  with  experiment  for  several  configu¬ 
rations  indicates  that  accuracies  of  ninety  percent  or  better  can 
be  obtained  for  force  coefficients  of  most  configurations.  This 
is  at  a  cost  of  about  $30.  for  ten  Mach  numbers  in  the  range 

0  0^3. 

3.  A  second  order  axial  perturbation  solution  can  be  combined  with 
modified  Newtonian  theory  to  adequately  predict  pressures  on 
general  shaped  bodies  of  revolution.  This  is  true  for  Mach  numbers 
as  low  as  1.2  even  though  Newtonian  theory  was  derived  for  high 
Mach  number  flow. 

4.  A  first  order  inviscid  crossflow  solution  is  not  sufficient  to 
predict  afterbody  or  boattail  lift  at  low  supersonic  Mach  numbers. 
However,  when  account  is  made  for  the  boundary  layer,  markedly 
improved  results  for  boattail  lift  was  obtained. 

5.  There  is  still  no  adequate  theory  available  in  transonic  flow 
which  is  computationally  fast  and  accurate  and  can  consider  blunt 
nosed  configurations  with  discontinuities  along  the  ogive.  Thus 
more  research  needs  to  be  directed  along  these  lines. 
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FIGURE  1  TYPICAL  BOCY  GEOMETRY 
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FIGURE 2  BOUNDARIES  OF  PERTURBATION  AND 
NEWTONIAN  THEORY 
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FIGURE  3.  TRANSONIC  WAVE  DRAG  OF  TANGENT  OGIVES. 
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FIGURE  12  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  BLUNTED  CONE; 

rn/rB  =0.35IM0o=l.5,a  =  8° 


\\*v — 

Gk\  A  A  A 


(ft(DEG) 

0 


EXP.(REF.8)  THEORY 


SHOCK  EXPANSION  THEORY 


FIGURE  13  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  BLUNTED 

CONE;  rn/rB  =0.35,  Moo  *2.96,  a-6* 


% ERROR  ( 


FIGURE  1 4  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR 
A  BLUNTED  CONE;  M«oF  l.5,0c-IO». 


% ERROR ( 


FIGURE  15  COMPARISON  OF  THEORY  AND  EXPERIMENT 
FOR  BLUNT  CONE;0C  «  10°,  =0.2. 


XCp(CALS.  FROM  BASE) 


FIGURE  16.  COMPARISON  OF  THEORY  AND  EXPERIMENT 
FOR  BLUNTED  CONE;  ^  -10*  rp/r#«0.4. 
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Xcp/* 


j 


ANGLE  OF  ATTACK  (DEGREES) 

FIGURE  17  COMPARISON  OF  THEORY  WITH  EXPERIMENT 
FOR  TANGENT  OGIVE- CYLINDER. 

.1-14  CALIBERS 
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FIGURE  18  COMPARISON  OF  THEORY  AND  EXPERIMENT 
FOR  CONES  OF  VARIOUS  LENGTHS. 


8 


XcptCALS.  FROM  NOSE) 


FIGURE  1 9  COMPARISON  OF  PRESENT  THEORY  WITH 
EXPERIMENT  AS  A  FUNCTION  OF  AFTERBODY  LG. 
(2.83  CALIBER  TANGENT  OGIVE  NOSE). 
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FIGURE  20.  ZERO  UFT  DRAG 


THEORY 


I 


FIGURE  21  COMPARISON  THEORY  AND  TEST  DATA 
FOR  5m/54  RAP  PROJECTILE. 


CMa(about  C.G.) 


'Nja=0°) 


0  BRL  BALLISTIC  RANGE 
®  AEDC  WWD  TUNNEL 
a  NOL  WIND  TUNNEL 
—  THEORY 


MACH  NUMBER-Mop 

RGURE22  COMPARISON  OF  THEORY  AND  TEST  DATA  FOR 
IMPROVED  5^54  PROJECTILE. 


o  os  i1o  1*5 2*5  2T5 

MACH  NUMBER  - 

FIGURE 23  COMPARISON  OF  THEORY  AND  TEST  DATA  FOR  I75MM 

XM437  PROJECTILE 
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Xcp(CALS.  FROM  BASE) 


FIGURE  24.  COMPARISON  OF  THEORY  AND  TEST  DATA 

FOR  155 MM  PROJECTILE 


5 


Xcp( CALS.  FROM  BASE) 


FIGURE  25  AERODYNAMICS  OF  5  INCH  GUIDED 

PROJECTILE  BODY 
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GLOSSARY 


ca 

Caba 

CAf 

Carb 

CA 

Avis 

caw 

% 


c 


M 


°Ma 

CN 


d 


dB 

H 

i 


M 


Pr 

R 


R' 


Total  axial  force  coefficient 

Axial  force  coefficient  contribution  from  base  pressure 

Axial  force  coefficient  contribution  from  skin  friction 

Axial  force  coefficient  contribution  from  rotating  band 

Axial  force  coefficient  contribution  from  viscous  separation 
on  nose  at  subsonic  Mach  numbers 

Axial  force  coefficient  contribution  from  expansion  and 
shock  waves 

Zero  lift  drag  coefficient;  =  C,^ 

Crossflow  drag  coefficient 

Mean  skin  friction  coefficient  based  on  freestream  Reynolds 
number 

Pitching  moment  coefficient  about  nose  unless  otherwise 
specified  (positive  nose-up) 

Pitching  moment  coefficient  derivative  -  dC^/dct 
Normal  force  coefficient 

Normal  force  coefficient  derivative  -  dC^/da 
Pressure  coefficient;  Cp  =  (P-P^/l/Zp^V^,2 
Diameter  (calibers) 

Base  diameter 

Mean  height  of  rotating  band  in  calibers 
Body  Length  (calibers) 

Mach  number 

Prandtl  number 

Body  Radius  (calibers) 

dR/dx 
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R\[  Reynolds  number  -  (pV£)/iJ 

Rj  Turbulent  boundary  layer  recovery  factor 

Sw  Wetted  surface  area  of  body 

Sp  Planform  area  of  body 

Tw  Wall  temperature 

u,v,w  Velocity  components  in  cylindrical  coordinate  system 

V  Total  velocity  -  V  =  /u2+v2+w2 

Vol  Volume  of  body 

x,r,0  Cylindrical  coordinates  with  x  along  axis  of  symmetry  and 
in  calibers 

x,y,z  Rectangular  coordinates  with  x  along  axis  of  symmetry  and 
in  calibers 

xCp  Center  of  pressure  in  calibers  from  nose  unless  otherwise 
specified 

Xp  Distance  to  centroid  of  planform  area  in  calibers  from  nose 

xu,ru  Coordinates  of  point  below  which  perturbation  theory  cannot 

be  applied 

Xj  Distance  measured  relative  to  shoulder  of  boattail 

a  Angle  of  Attack 

8  Angle  between  tangent  to  body  surface  and  axis  of  symmetry 

Y  Ratio  of  specific  heats  (y  =  1.4) 

6  Angle  between  a  tangent  to  the  body  surface  and  freestream 

direction 

6*  Angle  which  the  nose  makes  with  the  shoulder  of  the  body  (degrees) 

£  Velocity  potential  in  cross  flow  direction 

T1  Ratio  of  drag  coefficient  of  a  circular  cylinder  of  finite 

length  to  that  of  a  circular  cylinder  of  infinite  length 

0  Cylindrical  coordinate  measured  with  9  =  0  in  leeward  plane 
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0C  Cone  half  angle 

y  Coefficient  of  absolute  viscosity 

p  Density 

<j>  Total  velocity  potential  which  is  made  up  of  axial  and 

crossflow  velocity  potentials 

Velocity  potential  for  axial  flow 
Subscripts 

00  Freestream  conditions 

a  Afterbody 

B  Boattail 

BA  Base 

n  Nose 

o  Stagnation 

r  Reference  conditions  (reference  length  is  the  afterbody 

diameter  =  dr) 
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COMPUTER  PROGRAM  TO  DETERMINE  PRESSURE  DISTRIBUTIONS  AND  FORCES  ON 
UNGUIDED  PROJECTILES  OR  THE  BODY  ALONE  OF  TIIF.  GUIDE!)  PROJECTILE 


The  methods  described  in  the  report  to  obtain  surface  pressures 
and  force  coefficients  have  been  programmed  for  high-speed  digital 
computation.  The  purpose  of  this  appendix  is  to  provide  a  general 
description  of  the  program  including  a  listing  of  the  program  and 
a  sample  of  the  required  input  and  resulting  output. 

A.  DESCRIPTION  OF  PROGRAM 

The  program  reads  in  the  body  geometry  with  x  =  0  as  shown  in 
Figure  1.  If  the  nose  is  truncated,  a  conical  nose  of  angle  given 
by  Figure  2  is  automatically  placed  on  the  truncated  portion  to 
get  the  perturbation  solution  started;  but  the  pressure  integration 
begins  at  x  =  0,  which  is  the  location  of  the  first  point  read  in. 

If  the  nose  has  a  spherical  cap,  then  the  program  automatically 
computes  this  and  again  the  first  point  read  into  the  computer  is 
at  x  =  0.  However,  the  pressure  integration  begins  at  x  =  -rn. 

It  is  suggested  that  the  description  of  the  body  be  read  in  to 
at  least  three  decimal  places  if  possible  because  the  resulting  solution 
will  not  be  as  accurate  as  it  could  be  otherwise.  For  example,  if 
the  ogive  has  a  formula  it  is  suggested  a  desk  calculator  be  used 
to  compute  the  body  coordinates  as  opposed  to  a  slide  rule. 

Once  the  coordinates  of  the  body  are  read  in  (the  various  body 
geometry  options  are  discussed  below) ,  the  program  then  computes 
a  new  set  of  body  coordinates  where  the  flow  field  solution  will 
actually  be  found.  These  points  are  unequally  spaced  along  the 
body  to  conserve  computational  time  but  are  also  soaced  closely 
enough  so  an  accurate  solution  can  be  assured.  Once  the  body 
geometry  has  been  found,  the  program  checks  to  see  whether  the  Mach 
number  is  subsonic  (Moo  £  0.8),  transonic  (0.8  <  Moo  <  1.2),  or 
supersonic  (NLo  >  1.2)  and  then  proceeds  to  numerically  calculate 
the  force  coefficients  for  that  particular  Mach  number. 

B.  INPUT  DATA  CARDS 


CARD  NUMBER 


PARAMETERS  READ 


FORMAT 


1 


M 


13 


2 

3 

4 

5 


AL,  DIA,  HB,  AINF,  RHOINF,  AMUINF,  (4F 10 . 4 , 2F15 . 12 , 
I PRINT  15) 


MN  (MN  <  16) 


13 


AM  (I) 

I  =  1,2,3, . . . ,  MN-1 ,  MN 


16F5.3 


N,  NSHAPE,  Nl,  N2 ,  N3,  NBLUNT  NFL,  (815,  4F10.5) 
NN1A,  C2 ,  C4 ,  F,  RR 
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CARD  NUMBER 


PARAMETERS  READ 


FORMAT 


6 


X(I),  R C I ) 

I  =  1,2,3,...,  N-l,  N 


2F15.10 


N+5 


C.  DEFINITION  OF  PARAMETERS 


PARAMETER 

M 


AL 
DIA 
HB 
AINF 
R1IOINF 
AMU INF 
IPRINT 


USE 

Specifies  number  of  cases  to  be  run.  If  M  >  1,  then 
only  one  data  card  needs  to  be  included  for  M  but 
cards  2  through  (N+5)  are  included  for  each  additional 
case 

Angle  of  attack  (degrees) 

Body  reference  diameter  (feet) 

Mean  height  of  rotating  band  in  calibers 

Freestream  speed  of  sound  (ft/sec) 

Freestream  density  (slugs/ft3) 

Freestream  absolute  viscosity  (lb-sec/ft2) 

IPRINT  =  1;  pressure  coefficients  are  to  be  printed 
2;  no  pressure  coefficients  printed 


MN 


Number  of  Mach  numbers  where  solutions  are  computed 
(MN  <  16) 


AM  ( I ) 


Mach  number  where  solution  is  computed 


N  Total  number  of  points  read  in  along  body  surface 

(N  <  30) 
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PARAMETER 


USE 


NSHAPE 


N1 


N2 


N3 


N BLUNT 

NFL 

NN1A 


Parameter  which  describes  body  shape 
1.  Pointed  Body 


NSIIAPE  =  1;  nose  only 
2;  nose  plus 
3;  nose  with 
There  may 
afterbody 
4;  nose  plus 
5;  nose  with 
afterbody 


afterbody 

discontinuity  in  it. 
or  may  not  be  an 
present . 

afterbody  plus  boattail 
discontinuity  in  it  plus 
plus  boattail 


2.  Blunted  or  Truncated  Nose 


NSHAPE  =  3;  nose  with  or  without  discontinuity. 

There  may  or  may  not  be  an  after¬ 
body  present. 

5;  same  as  above  except  afterbody  and 
boattail  are  present. 

If  NSHAPE  =  3  or  5,  at  least  5  points  must  be  read 
in  along  each  of  the  ogives  even  if  the  ogive  is  a 
straight  line. 

Number  of  points  read  in  along  first  ogive 

Number  of  points  read  in  up  through  the  second 
ogive  (includes  first  ogive) 

N3  =  1;  conical  boattail 

2;  ogival  boattail  (at  least  5  points  must 
be  given  along  boattail  if  N3  =  2) 

NBLUNT  =  1;  pointed  body 

2;  truncated  or  spherical  cap 

NFL  =  1;  spherical  cap 
2;  truncated  nose 


NN1A  =  1;  Blunted  nose  with  no  discontinuities 
present  other  than  the  intersection 
of  the  nose  cap  with  the  ogive  (N1  =  1 
and  N2  >_  5)  . 

2;  Blunted  nose  with  a  discontinuity  in 
the  ogive  so  there  are  two  ogives 
present  (N 1  >  5  and  N2  >  9). 
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PARAMETER 


USE 


C2,C4 


F 


RR 

X(I),R(I) 


Parameters  which  specify  mesh  spacing.  If  nose 
is  pointed, C2  =  0.9  and  C2  =  20  are  nominal  values. 
For  other  nose  shapes,  C2  =  .05  and  C4  =  1.0  are 
nominal  values. 

Constant  which  determines  limiting  body  slope  for 
a  given  Mach  number.  F  <_  1.0  with  F  *  0.95  recom¬ 
mended. 

Radius  of  spherical  cap  or  truncated  meplat  in 
calibers . 

Body  coordinates  (in  calibers)  where  I  <  30. 


D.  PROGRAM  LISTING 

The  Fortran  listing  of  the  source  desk  currently  being  used  at 
the  Naval  Weapons  Laboratory  is  as  follows: 
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1 


fft 

I 

U1 


PROGRAM  “AIN  TRACE 


CrC  6600  FTN  V3.0-HA06  OPT=0  u9/12/’2  If. 20. 07. 


5 


10 


15 


20 


25 


30 


35 


GO 


45 


*0 


r  5 


PROGRAM  MAIN(CUTF'JT  ,lNPUT,TAPEf  =  I NPU T  , T A PE6=0UTPUT ) 
COMMON/GEOm/RP(6),XC30>,R<30) ,C2,N.NSHAPE,N1,N2, XO  ( 225) , R3 ( 2 25 > 
CO“MOM/G5C1/  REP  (225)  , DETA 

C0MMON/GEO2/NM  , NN 2, NN 3 , NN4 , NFL , NBL UNT , NN , NNI ,  IPRINT.NN1A 
C0MM0N/GE03/V0VS,AL,XM,YM,XINT,YINT,NMA 
COMMON/GE04/K, F ,90 ,  RREF 

C  CMMON/ 91 S2  /  SLM1,SUM2,SUM3,SUM4,SUM5,SUM6 

COMMON/ 9 ATI/  T (1 00 > ,AK< 100) ,  AE(inC) ,C (225) ,C1 (225> ,C3 

C0MM0N/9ISC/  I,JK,AI2,SUM, JH.PI 

COMMCN/9ASE/CAP,CN0,CM!3 

COMMON/ 3 ANC/CAF,CNP,CMP,H9 

COMMOM/OIS1/  J 1 , J3 

C0““0N/WAVE/CA9L,CNBL  ,  CMOL  ,  C  A  * ,  C.N  H,  CMW 
COMMON/ VOL/  VOL ,CAF ,  C NF , CMF , R N , 9 1  A, XP, 4®, VOLN 
COHMON/ICOU/  ICOONT 

DIMENSION  AM  (20) ,CN<20> ,CM(20) ,CL(20)  ,CO(20> .XCPI20) ,CNAL (20) , 
1CMAL (2  0 ) ,CA1 (20) ,CAF1 (20) ,C4B1(20) ,CAW1(20),C4P1 (20) ,ETA(9> , 
2AL0019) . AMC (IP) ,COC (1 0) 

OATA(ET4(I),I-l,7)/.s3».5Ti.613j.6G..E65».70>.765/ 

DATA(4LOO(I),I=l,7)/l.,2.,t.,6.,B.,12.,20./ 

OATA(AMC( I) ,1=1,9) /0. ,.3, .4,. 5,. 7,. 9,. 9, 1. ,1.4/ 

0 ATAICDC (I) ,1=1,9) /l. 2, 1.2, 1.25, 1.35, 1.74,1. 82, 1.82, 1.8, 1.53/ 

READ ( 5,50)  M 
50  FORMATII3) 

C  M=NUM9ER  OF  CASES  TO  BE  COMPUTED. 

90  27  MM=1,M 

RE  40(5,4  3)  4L,0IA,HP,AINF,PH0INF,  4MUKF  ,  IPRINT 
43  FORMAT (4F10.4,2F15.12,I5> 

C  AL=ANGLE  OF  ATTACK  (DEG)  D I A  =  REFE  R  E  NCE  DIAMETER  OFF  (JODY(FT). 

C  AINF.RMOINF.AMUINF  ArE  ThF  FPEcSTREAM  REFERENCE  CONDITIONS  F  Ok 

C  SPEED  OF  SOUND (FT/SEC)  ,  DE NS  I T Y ( SLUGS/F T ” 3 ) , AND  ABSOLUTE 

VISCOSITY (LB-SEC/FT”2)  RESPECTIVELY  AT  THE  GIVEN  ALTITUDE 
C  IPRINT=1  IF  PRESSURE  COEFFICIFNTS  ARE  TO  3E  PRINTED  -2  OTHlRWISE 
C  H«=ME4N  HEIGHT  OF  ROTATING  BANC  IN  CAlI3ERS,  IF  NO  RANO  3JESENT  HB  =  9. 
WRITE (F, 6)  mm,AL,OIA 

6  FORMAT!//, FOX, ’CASE  NO. * , 13 , //,3 0 X , • ANGLF  OF  ATTACK  =’,Fb.2, 

1*  QEGS’ ,10Xf*REF£PENCc  DIAMETER  = * , F6. 3, *FT*, //) 

WRITE (6,7)  AINF.RMOINF.AMUINF 

7  FORMAT (54X, ’REFERENCE  CONDIT IONS*  , //  ,54X, ’SPEED  OF  SOUND  =* , 

1F9.3,’  FT/R£C*,/,54X, ’DENSITY  =*,F10.7,*  SLJGS/FT”3 

2  *,/ ,54X,*APS0LUTE  VISCOSITY  =’,F15.1?,*  L B-SEC/FT  ”2* , / / ) 

AL=Al/c7.295A7 

ICOUNT  =  0 

RE A DIE ,50 )  MN 

MN=NUM PE R  OF  MACH  NUMBERS  TO  COMPUTE  THE  FO^Cf  COEFFICIENTS  CF 
C  A  PARTICULAR  CASE. 

3EAO ( 5 , 1 5 )  (AM(I),I=1,MN) 

15  FORMAT (16F5. 3) 

DO  1  J=1 , MN 
ICOUNT=ICOUNT*l 
VOVS-  Am ( J ) 

PRFF=".5 

°ETA  =  SPRT  (ARS  I  VOVS”2-l.  )  ) 

IF ( BET  A.LE.0.5)  RFTArO.5 
CALL  Gr 9 M 
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PPOGPA" 


60 


65 


70 


75 


80 


35 


90 


95 


100 


105 


MAIN  TRACE  CDC  6600  FT  N  Vi. 0-0306  OPT  =  Q 

A  M ( J)  =  VO VS 
IF  CJ.GT.l)  GO  TO  17 
10  IFCN1.NE.21  GO  TC  17 
THEC=ATAN(RP(1)) 

THETA=THEC*57 .29583 
WRITE (6 | 301  TMETA 

30  FORMAT (IX, 17HC0NF  HALF  ANGLE  =,F10.5,/) 

17  CONTINUE 

VINF=VOVS*AINF 
9  N  =  PHO IN  F* V  I NF / AMU  I NF 
CALL  SKINF 
CALL  8ASEP 
CALL  9BAND 

IFCAL.LT. 0.0001)  GO  TO  18 
IF(VOVS.LT.1.2>  CALL  NORMFO 

18  IF(VOVS.GE.O.til)  GO  TO  19 
ICT=NN1 

IFINSHAPE.E0.3)  IC  T  =  NN2 
IFCNSHAPE.E0.5)  ICT=NN2 
THE1=ATAN(R8P (ICT) ) ’57.293 
IFCTHE1.GE.10.)  GO  TO  51 
CAW=0 . 

GO  TO  c 

51  CAW=0.012*(THE1-10.) 

GO  TO  5 

19  IFCVOVS.LT. 1.19)  GO  TO  2 
CALL  HYBRID 

GO  TO  5 
2  CALL  TRANS 
5  CA=CAF*CA9+CAN+CAP 
CA1CJ)=CA 
CAF1 ( J)=CAF 
CA91CJ)=CA8 
CAH1C J)=CAH 
CAP1(J)=CAP 
XT  =XB ( NN  > +RR 

CALL  INTERD ( ALOO , ETA , XT , ET A1 , 7 , 3 ) 

APEF=3.14159*RREF**2 

AMCl=VOVS*SINCAL) 

CALL  INTERP (AMC.COC.AMCl ,C0C1,9, 3) 

CNV=CDC1»ETA1»AP*AL’*2/AREF 

CMV=-ETA1*CDC1*AF*Al**2*XP/(APFf*2.*RREF) 

IFtAL.GT. 0.0175)  GO  TO  52 
CNV=0. 

CMV=0. 

52  CN(J)=CNF+CNB+CNk*CNP*CNV 
CM ( J)=CMF+CMB*CMW*CMP+CMV 

CL (J>=CN(J)*COSCAL)-CA»SIN(AL) 

CD(J)=CNCJ>*SINCAL)*CA»CCS(AL) 

IF (ABSCALJ.LT.O. CO 01)  GO  TO  1 
CMAL(J)=CMCJ) /AL 
XCPCJ) =-CM(J> /CN  (J) 

CNALCJ)=CN(J)/AL 
1  CONTINUE 
WRITr (6,8) 


3®/ 12/72  16.20. D7. 
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/ 


4 


110 


Ovor-P'M' 


115 


120 


125 


MflTN  T,4rC  l^T  *>f>00  FTn  •!  % .  u  — ^30 «  OPT-u  u«/lc/72  lb. 20. 07. 

«  F  ^“AT(//,53X,«flnflL  F0»CE  COF'TK!nUirCKS*,//,lX,»iaCH  N0.*,lfcX,*SK 
iin  fp  iction*  ,  i4x  ,*t’ASE  pressu^i* , i xx, * pressure* , i4x,*ppct rusicns* 

2»1<*Xi* TOTAL*,//) 
no  ’1  L-  1  ,MN 

WRITE <6, 9)  A'*<L),CAFlIL),CA91(L),€AWl(.),CAPl(t),CAl(L) 

9  F0RMAT(3X,F6.3,18X,F6.4,2'JX,F6.<.,17X,F6.4I16X,F6..,,17X,F6.4) 

31  CONTINUE 

WRITE (6,12) 

12  FORMS T(// >E6X, "FORCE  COEFFICIENTS*, //,10X,*MACH  NO . * , 1  OX , *CD* , 

11CX,*CN",10X,*CL*,10X,*CF.  *,10X,*CNAL*,10X,*CMAL*,10X,*XCP/D*,//) 
no  14  L  =1 , MN 

WRITE  (6, 13)  AM(L),CD(U,CN(U,CL(L)  ,CM(L)  ,CNAL  1  L)  ,CMAL(L)  ,  XCF  (U 
14  CONTINUE 

13  F‘0RMST(12X,F5.3,9X,F6.4,6X,F6.4,tX,F6.4,6X,F6.3,6X,F6.3,8X,F7.3, 

10X.F7.4) 

27  CONTINUE 
E  NO 
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FUNCTION  «R'FCH<7) 

A5>SECH=4LCG  (1  ./Z*SCRT  (1  ./7*«  2- 

5FTURN 

CNO 


SU'TRPUT  I N:  BflSco  TRACT 


CJC  bCHT  FTN  v3.n-»30e  Of  T=  0  ;o/i 2/7c  lb.2C.37. 


5 


10 


15 


?0 


S  U3F'3u  T INC  rtSTP 

;cmm0n/'G!’OM/RP  (6)  ,*<30> ,R<3D)  ,C2>!l,MShAPF,Nl,N2,XR(225>,K9(22E) 
COM^M/GECi/  RRF<??5>  ,BETA 

■■•0MM0N/GC02/NN1  ,NN2  ,NN3,  NNR,  NEl,  NTi_UNT  ,NN,  NNl ,  1  PRINT, NN1  £ 

COMMOM/r,EC’/VCVS,AL,XM,YM,  XINT  ,  YTNT  ,NMa 

COMMON/GEOW  K,F,PR,RREF 

COMMON/TSSP/  CA3,CNS,CM9 

DIMENSION  TTP3O<20> , TXM1 (20) 

(]ATA(  TCP 30<I)  ,  T=  1,19)/.  120,.  12*.,.  13  0,.  135,.  1A2,.  154,.  209,.  219,. 221 

I, . 218, .211, .191, .173, .157, .143, .131, .11*,, .104, .095/ 

0ATA(TXM1(I)  , 1=1, 19>/0.,. 5, .7, .8, .8=, .90, 1.0, 1.05, 1.1, 1.2, 1.3, 1.5, 

II. 7,1.9,2.1,2.3,2.5,2.8,3.0/ 

Cfl  4W=  0  . 

THIS  SUBROUTINE  CALCULATES  BASE  DRAG  THROUGHOUT  THE  MACH  NUMBER  RANGt 
0CPP£=0. 

cau  INTTRP(TXM1,TCP3D,V0VS,CP30,19,3) 

CDRP  =  CP30* (RPCNN1/RREF) **3 
IFTat.LE. 0.0175)  GO  rc  3 

0CPoa=(.312-.0036*0OVS)‘aL*57.295»(9B(NN)/RREF)»»3 
3  CONTINUE 

cae=  CD8P  ♦  ocpea  .  caaw 
c  n  n = o . 

CM9=0 . 

RETURN 

END 


°AGE  1 


25 


01-8 


SUiSPOUTINE  BLUM  TRACr 


CPC  *>600  FIN  U3 .0-Po0e  OPTO  P9/12/7 i 


SUBROUTINE  SLUM 

GCHMON/GEOH/PF  (F.)  ,X(33)  ,P(?Q)  ,C2,N,NSHAPfc,Nl,N2,  XB(225)  ,RB  (225) 
GOHMON/GEOl/  R8P(225> , BETA 

COMPCN/GFO?/NM,KN2t’IN7,NNl.,HFL,NbLUNT,NN,.NNI,IPRINT|NNlfl 
COHMON/GF03/VCVS,AL,XM, YH , XI  NT , Y I  NT , NNI  a 
COMMON/GEOA/K.F.CR.RREF 
0 OMMO  N/GE  05 /  C 3 

10R5)  F°5  <X<1>  »X(1)  ,X(2)  ,X(3)  ,X(4)  t  X  (  5  ) >  R  ( 1  )  ,R  (2)  ,R  (  3)  ,  R  (  4  )  ,  R  (5  )  , 

10  VOV=VOVS 

IF (VOV'.LE. 1. 19)  VOV  =  l.O0()l 
AMU=ASIN(1. /VOV  )*F 
IF ( ORB • L  E  •  A mU )  GO  TO  21 
A  MU=F  »DP9 

15  VOVS=l./SIN (AMU) 

BFTA=SORT (ABS(V0VS*»2-1.) ) 


21 

J  =  1 

TH1  =  AT A*J  (ORB) 

RR-R 1 1 ) »COS (THl) 

20 

2 

TFIVOVS.GT.2.1)  GO  TO  14 
THET1=27. 5/57. 295 

OMAN  (THETl) 

IF(D.LT.DRB)  0=DRB 
THFtl=ATANIO) 

25 

XM--RR»SIN(THET1) 

YH  =  RR’COS (THETl) 

XB (1) =XM-YM/TAN(THETl) 

GO  TO  15 

1 A 

YH=RR*BETA/VOVS 

30 

XHs-RR/VOVS 

XB(1)=PR»»2/XM 

THE  =ATAN(YH/(XM-XB(1)) ) 
THETl  =  TH  E*  F 

XM=-PP*SIN( THETl ) 

35 

YH=  PR*C  OS ( THET 1 ) 
XB(1)=XM-Y"/TAN( THETl) 

15 

RP(1) =TAN(THETl) 

T  HEM  HE  Tl*57.295 

R9  (1) =0. 

A  0 

XB(2)=XM 

RB(2)=YH 

R9P(1)=TAN(TmET1) 
R9P(2)=RRP(1 ) 

Z=SQRT (1. *0Rp**2) 

45 

XI=-OR?»RR/Z 

R I=RR/Z 

X INT  =  X  I 

YINTrRI 

XIXH=APS(XI-XM> 

50 

HNT  A  =  2 

IF(XIXH-.OCl)  ie,18,lP 

15 

K  =  1 

N  N 1  =  2 

GO  TO  15 

55 

19 

F=15./V0V,»? 

lb. 20. 37. 


°AGE  1 


SUBROUTINE  ->lUNT 


t-  ■irr 


r^r  fcn-jQ  rr.j  vj.o-Pjj;#  a i- t « •: 


D9/'c/72  16.2J.07 


RAGE 


60 


65 


70 


75 


“0 


85 


go 


95 


ISO 


1  85 


no 


l-'iE.u.n  e=3. 

XT  (7)  =  »  ]  (  J)  .  ,  r  1  /  VO  V  *  *r*  Ru  (2  )  *  (Tr'CV-i '  . )  *  *2 
rn  7  <=j,i80 

a -k-p 

3. 9(0=  SORT  (Po^j-XPdO  ’•  21 

XB(K+1)=XG(KJ  ‘.Ol/VOV  »*r»°a<K) •a«*C.5U* (THE/30. )**2 
5.0P  (KJ=-X"(K)  /=-0(K> 

IF(X9(K*1J.GE.XI)  r.C  TO  10 
3  CONTINUE 
IQ  X9(K»1)=XI 
R”  (<*1  >=R1 
R9P(K*1> =rR9 
IF (NN1A, EC. 2)  GO  TO  16 
NN1  =X  » 1 
NN2=NN1* 1 0 
K=K  +  1 

XB(K*1)*XI 
R8  (K*l )  =  RI 
R9P(K«1)=9PP(K) 

IF(NFL.E0.2>  NNI=K*1 
IF  (RI.GE.QR)  GO  TO  99 
16  0X=-XI/6. 

09= (R (1) -PI) / 6 . 

00  13  J= 1 , 5 
K  =  K  +  1 

X3(K*1)=X3(K)+DX 
R9(K‘1)=R°(K)»0R 
ROP ( K  * 1 ) =  0R9 
13  CONTINUE 

I F (NFL . EQ  .  ?  )  NNI =K ♦ 1 
IF (NFL.EQ.2)  NN1=<*1 
IF(NNia.FQ.l)  GO  TO  99 
K  =  K*1 
IJ=1 

X9(K*l)=X(l) 

RB(<+1)=R<1> 

oact  FC°6 (X ,P, xe (<+l ) tRPP(K*l> ,N1 ,3) 

20  9ETl  =  9'T,l 

IFI9ET1 ,9T.  1.)  6eTl  =  l. 

<  =  <♦! 

XP(F+l)=X9(K)*C3*RETl*Ra(KJ 

CAU.  INTE»F  (X,R,X3(K+1>  ,RO(Ktl)  ,m,  3) 

CALL  FD»5(x,-,X!MK*l>  f  R-3D  (K+  1 )  ,*J  1 , 3 ) 

17  <  =  K  +  1 

IJ=IJ+1 
A  =  I J 
C5=A*C2 
0ET:='’=Ta 

I- (0ET1.GT. 1. )  9ET 1  =  1  . 
xD(K*!)=X9(K) tC5»BFTl*RB(K) 
iF(x9(xn).Gr.x(M))  <n(x*ii«x(Nn 
CALL  INTERF  (X,.7,X9  (K»l)  ,R  1  (<*1)  ,M  ,  3) 

CALL  R'PP  (V»PjX:7(K  +  l  >,R-JP(<«1)  ,f.'!,3) 
rc  1 »"  (-•  j . lt  .  (X(mi  , -. 03G 1 1 1  *.(>  rr  17 
[  J=  T 


2 


SJTPCUT  IKF  'M-UNT  T?3Cf  COC  eb03  FT  N  V5.0-P303  U<=1  =  0  li^/lZ/Td  It. 23.07,  PAGE 


•S'.' 


c 


10 


10 


:>c 


cf  Jj  "Tfi  «,HMs  0PT-3  00/12/?2  It-. 20. 37. 

•i'jppr;uTr>i-  crsci 

COM«0'I/.V  r'H/tr-(eif>  (3  C)  ,  '  J  3) ,  ?,N,  \5“A  =  '-,Nl,M2,X'l(22«i)  ,R1  (t\U> 

“OMHON/GCri/  ?DP (2?c)  ,  p  E  T  A 

CO-HCN/GFOLV'  N  1  ,NM2,-JN3,t.-.u,nrL,  (.oLUMT  ,  NH,  NNT  ,  I  POT  NT  ,  \'N1  A 

COM-fCM/ITPr/I  ,  J<,ai2,?ui,  Jr.,  PI 

CO“MCN/3A  T1  /  T  (ICO) ,AK(1'0) , AE (100)  ,C (22  =  ) , Cl  (225) ,C3 
0  =  5'>S  (POP  (JH)  -»B°(  JH-1)  ) 

XI  =  V0(JH)-!>ETfi*!?n(JH) 

T flU= BETA*  (!>/•(  X ^  ( I  >  -XI) 

IP  (TAU.3E.  l .  )  TAU=0.AiS'3'39')t> 

CALL  INT‘3y°<TfAK,TAU,AKX,10  0,3) 

CALL  INTr~>D(T,AE,iAU)Airy,100,3) 

C  CUEVA^E  SOLUTION  FOR  FI»ST  OPDEP  FU)CTIQN 

THP2=C ( JH) *8ET  A* SORT (XB ( I ) -<I ) *A . /? .  *Ef33T (?. ) /PI »SQRT  ( 1,  +T  AU) ‘ 

1  (AtX/TAU-AKV)  *  SQ  ( 62  T  A  *  PI  (JH)  ) 

HALF=0  . 

IF  (D.LT  .0.2  PCD  GO  TO  06 

HALF=C  (  JH)*2.*0ETA/PI*iQOT  (OB  (JH  )/FQ  (I)  )  *  SO*  T  (  ?  .  *  T  A(J  /  ( 1 . 4  T  AU)  )  * 
1<(1.+TAU)/TAU»AEX-AKY) 

9«  SUH=SUf'*-THF?+HALE 
RETURN 
ENO 


PAGE  1 


33 

I 


Oo 


B-14 


^UTOOU’If^  CJ I  S'' 


TRirr 


lDT  ofcflO  FTn  V/3.0— “iOf  UFT  =  0  09/12/74  lf.2tl.07. 


10 


15 


?0 


?5 


30 


35 


1*0 


45 


99 


SUOFCUTIMf  Cisr.2 

C0M“CN/GE9M/-  P (6) ,X  <  301  ,9(70) ,C2,M,NSfSPf,Nl,N2,Xy(22,),Kn(?2fc) 
COMMON/GF 01 /  PBP(225)  ,  r>ET  A 

Cn«M0N/94Tl/  T  (100)  ,  A  K  (  lu(i)  ,  A  E  ( 1 L  0 )  j'"  (22  51  ,Cl  (225)  ,C3 
CC9MCM/9ISC/I ,JK,fll?,SUM, JH,  PI 
C09M0N/9IS7/  SG“1,SU92,SU"3,SU“4,SU“5,SUH6 
CURVATURE  SOLUTION. 

XI=XR(JH)-9ETA*RS(JH) 

T  AU=Br TA»  RB (I)/(X0(I)-XI) 

IP ( TflU .5E .1. 1  TAU=0 .999999° 

CALL  I  NT  FRP (T , AK , T  AU , AKX, 100 ,3) 

CALL  INTERF(T,AE,TAU, AEX, 100,3) 

A = SORT (X9(I)-XI) 

B=SORT ( 1 . ♦ T  A  U 1 
n=?./PI‘SORT (2.  ) 

A1=C(JH) ‘SORT (PETA*P9(JH) 1 

31 =SQRT (  ( 1 , ♦ T  A  U 1 *RB ( I ) / ( T  A  U*  R0 ( JH 1 ) ) 

SU91=SU91-A1  •S**1.5*4./9. *0*8* l ( 3 . ♦ I AU) » AKX-4 . » AEX) 

SUM2=SU9?-Al  *A*2.*D*3*  (AKX-AEX) 

SUM3=SUM3  *A 1  *3ETA*A*2./j.«o»8*(AFX/TAU-AKX> 

SUM4r  SU94 -A  1  /A*0»AKX/B 

SUM5=SU° 5  *A  1  *EETA/A*0/B* (  (l.  +  TAU) /TAU*AEX-AKX) 

SUM6=SU96-A1  *9ETA**2/A*(VB*(2.*(1.*TAU) /T AU** 2* AEX- ( 2  .  - T AU ) / 

ITAU’AKX) /3. 

E=ABS(R9P(JH) -RBF(JH-l) ) 

IP (E.LT. 0.0001)  GO  TO  99 
CORNER  SOLUTION 
A2=A**2 

P=SQRT(R9(JH)/RB(I) )/PI 
G=SQRT(2.*TAU>/° 

IFdAU.GT.  0.999)  GO  TO  2 
H=l./ (l.-TAU) 

«1=C( JH) 

SUH1  =  SUM1-H1*4.*A2*F*B**2»G»(AKX-AEX> 

SUR2  =  SUM2-H1»2.*F*G»A<X 

SUM3  =  SUH3*H1*2.*BET  A • F  *  G* (3**2*AEX/TAU-AKX) 

IF (TAU.LT .0.999)  GO  TO  1 
SUM4  =  SUM4*-H1/ (8. *BFTA*RR (JH) ) 

SU95=SU95*3.'hl/(8.*R9(JH) ) 

S(JM6=SU96-7.»BFTA»H1/(».»RF(JH)> 

GO  TO 

SUH4  =  SUM4*Hl*F/A2*h*G*(AKX-AEX) 

SU“5  =  SUM5*H1*BETA*F/A2»H*G*(AEX/TAU-A<X) 

SUM6  =  SUM6-H 1 *  BET  A  * *  2 /A 2* F * H* G*  ( < 2 . - 1 AU* *  2 ) /TA U* » 2* AcX- ( 2 . - T AU) 
1 /T  AU*  AKX ) 

RETURN 

FNO 


PAGE  1 


10 


15 


20 


25 


-I:>-  r-:r.  6500  5TPi  *3.3-03118  OFUO  C9/lc/72  15.20.37. 

S'JW'UTINF  CI'-'C"? 

C0p“'C'-1/GE0m/-'F  (6>,X(30),RC*1)  ,T2  ,  N,  NSHAPF ,  N1  ,N2  ,  X3  (  225)  ,  P<?(  i  25 ) 

COmmON/GEC1/  ">!!=  (2?S)  ,  BETA 

CO-*“CN/OAT17  T(iro),SK(UO),fl£(lCO)  ,C(225)  ,C1  (225»,C3 
C0MM0N/3ISC/I  ,  JK, A 12,  SUM,  JH.PI 
0=49S (R3P ( JH)  -RBC  (  JH-1)  ) 

XI=X!3  ( JH)  -?FTA*R3  (  JH) 

TaU=HFTJ*Rn(I)/(xq(I)-XI) 

IF  (TAU.GE.l.)  TAl=n. 9999999 
CALL  INTERP (T, AK.TAU, AKX.1C0 , 3) 

CALL  INT£RP(T,ftE,TAU,aEX,lC3,3) 

C  CURVATURE  SOLUTION  FOR  CCMPL IMENT ARY  FUNCTION 

A=SQPT(2.*TAU*R3(JH>/ (RB(I) * (l.+TAU) ) ) 

3  =  C 1 ( J H ) *  9E T A /c I ’SORT (  3E T 4 *R 0 ( JH > ) 

SUH^SUM  +  2 .*Bt TA*C1 ( JH )  /P I “ A  *  (  ( 1 . ♦ T A U ) /T AU* AEX- AKX ) 

IF(O.LE.O.COBl)  CO  TO  99 

C  CORNER  SOLUTION  FOR  COMPLIMENTARY  FUNCTION 

31=-C3*3ETA/PI 
IF (TAU.3T .0 .995)  GO  TO  1 

SUM1=9*SQRT (XP(I)-XI) *4,/3.*SQ=T(2.*(l.*TAU>)*(AEX/TftU-AKX> 

SUM2=31/ CXB (I) -XI)*A/ (l.-TAU)* (AFX/TAU-AKX) 

GO  TO  2 

1  SUM1=0. 

S  UM2= -3 . *  C3  /(8.*f9(I)  > 

2  SUM  =  SUM*-SUM14SUH2 

99  RETURN 

END 


RAGE  1 


SJ^OUTINF  i)ISC4  TRACr 


Cjr  6bCu  FTN  \/3.4-R30e  OPT  =  Q  G9/12/7c  16. 2d. 07. 


5 


10 


15 


20 


?5 


00 


05 


40 


C 


c 


1 


2 


SU'l°nUTINE  515-4 

COM«ON/GrOM/BB  (6) ,  X  (  5  0)  ,R  (  30)  , C2 , N, N5HAPE , N1 , N 2 ,  X9  (225) 
COUMON/GEC1/  PBR(225> ,3ETA 

C3MM0N/DflT 1 /  T  (100) ,  A  K  ( 1  n  9  ) ,AE(1CC) »C  l  2  25) ,C1(?25) ,  C  0 
COMMON/0 1 SC /I , JK,ai?,  SUM, JH.PI 
COMMON/D  IS2/SU"1,SUM2  ,SU“3  ,  SUM*- ,  SUM5  ,  SUM6 
CURVATURE  SOLUTION 
X  I  =  X9 (JH) -SET  A  *  R  5  (  JH) 

TAU=9FTA*Re(I)/(X3<I)-XI) 

IF (TAU. GE .1 . )  TAU=C. 995099° 

CALL  IMERP<T,AK, TAU,  AKX, 100,2) 

CALL  WFRF(T,AE,TAU,  AEX,1C0,3) 

A  =  SQR  T (XB(I)-XI) 

°=SOPT (1 . +T  A  U ) 

5  =  2.  **1. 5/PI 

E  =  AQS (R3P (JH) -R°P( JH-1)  ) 

F=SQRT(R3(JH)/R3(I))/PI 
G=SORT <2. *TAU) /8 
A  ?  =  A*  *  2 
H 1 =C1 ( JH) »F*G 

SUMl-SUMl-Hl*4.*A2*e*»2» (AKX-Ai X) 

SUM2=$UM2-H1'2.*AKX 

SUM3=SUM3*H1»2.«SETA» ( 6** 2* A EX/T A U- AKX ) 

IF  (F. LT.  0 .00  0  1  )  GO  TO  2 
CORNER  SOLUTION 

H2=C1 ( JH) *D*G»SQRT < BE TA*R9  < JH ) ) 
SUMl=SUHl-H2*A»»1.5*4./9.* ( ( 3 .4-T A U> * AKX-4. * AEX ) 
SUM2=SUM2-H2*A‘2.» (AKX-AEX) 

SUM3=SUH3»H2*eETA*A*2./3.* ( A fc X/T AU- AKX ) 

S  UM1  =  SUM  1  *-C  3  '2.*G*F»AKX 

IF (TAU.LT.O .995)  GO  TO  1 
SUM2=SUM2-C3  /(8.'8ETA*RB(JH) ) 

SUM3=SUM3-3.*C3  /<8.»R8(JH)) 

GO  TO  2 


H  =  1 . /  (l.-TAU) 

SUM2=SUM2-C3  /A2*F*H*G*  (AKX-AEX) 

SUM3=SUH  3 -C3  *BFT A/ A 2* F* H«G* < A FX/T AU- AKX) 

CONTINUE 

RETURN 

END 


, PR (225) 


RAGE  1 


Stn=C'JTIM5  FH5  t>ice 


c^r  J69C  STr.  1/.S.0-=>316  OFT*i)  Jb.20.3T. 


5 


la 


15 


S  U?pCin  I  ME  FC‘  (A,  X  1,  X?,X3,  Xi.  ,XE  ,  FI  ,F2  ,  f',  FU.FF,  FX) 


a  i  = ( x- >  ^ ) »(x-xc)*(2. 

a2=[v-X4> * ( X - X  5 ) *(2. 
a3=(x-x‘*)*(x-x5>*(2. 
ai.:(x-x5)»(x-xc)*(2. 

A5-<X-X3)  *(X-XM*(2. 
m=(Xl-X2)»  (X1-X31*  (Xl-Xi.)  *  (X1-X5) 
32=  (X2-X1>*  (X2-X3)  •  (X2-X1.)  * (X2-XE) 
D3=(X3-X1»* (X3-X2) *  <XJ-Xh) "(X3-X5) 
D4=(X1.-X1)»  (XI.-X2)  •  (X4-X3)  •  (XI.-X5) 
P5=(X5-X1 )  *  <  X  5 *X2 ) * (X5-XJ) » (X5-X4) 

ci=Ai/ni 

C2=A2/D2 

C3=A3/P3 


•X-X2-X3)*  (X-X2)*(>-X3)*<2.  *X-X-.-X5) 

‘X-X1-X3I  *  (X-X1)»(X-X3)’  (2.  *X-X*-X5> 

’X-Xl-X?)  KX-XU*  (X-X2J*  (2.  *X-Xh-X5) 

•X-X1-X2) ♦(X-X1J*(X-X?)*(2.  »X-X3-X5) 

•X-X1-X2) * (X-Xl >*  <X-X2>*  <2.  'X-X3-X4) 


C^sAA/Hi. 

C5=A?/P5 


FX=Cl*rnC2*F2+C3*F3*C4*F4*-C5*F5 

RETURN 

ENP 


0  1 


Q 

0 

Q  7 
0  9 

Q  9 
0  1U 
Q  11 
Q  12 
Q  13 
Q  1 
Q  15 
Q  16 
0  17 

a  ie 

a  19- 
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l^Hr 
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It. 2  0. 07 


PAGE 


SU^OUTINP  FTCSCTX.TY.X  .Y.N.J) 
niPFfiSION  T  >  (30)  ,TY  C30I 
I  =  n 

l  1=1+1 

5  IF  ( T  X  ( I )  ■  LE  .  X  >  OF.  TO  1 

K  =  J+2 

IF(I.LF.K)  I=K 
IFCI.GT.  ( N  ~  ? ) )  I=N-2 

CALL  c  05  <  X  ,TX(I-?)tTX(I-l»,Tx(I)fTX(r*l),TX(I*2»tTY(I-2)fTYlI-l> 
10  1>  ,TY(I),TY(I»1)  ,TYU*2>  ,Y  ) 

RETURN 

END 


33 

I 


O0 


vnpruTiN;  r.roM  t  »/■... 


•ib3C  CT»,  ^.o-pjoe  CF  T=  I  u“/l£/72  lfc.20.J7. 


p«GE  1 


s 


10 


15 


20 


tx> 


VO 


25 


30 


35 


1.0 


1.5 


50 


55 


su:,Po;jTiNr  g‘cm 

CO'*«OM/,ii  O“/»Fr6),X(30)  ,C2,N,NS)-4?F,N1,N2,'C9<225>,RB<22P> 

COMMcn/GEGI/  R  jF  (225)  ,3cTA 

"OHH0N/G5O2/WM  »NN2»NN3,NN‘.,NFL»NfcLtlNT  ,NN,  )NI ,  IP»iNT,  MM  ft 

C0MM0N/GEDW0VS,4L,XM,YM,XIN1  ,YINT  ,M  IA 

C0.MM0N/GE0<t/K  ,F,Ro,SRFf 

COMM0N/GS05/  C3 

00MM0N/IC0U/  ICOUNT 

C0mmCN7lE'IG/'>L,AM-  ,A_A 

IF  (ICOUNT.GT.l)  GO  TO  31 

9240(5,1)  N,NSHAFE,N1  ,N2f  N3,N3LUNT,NFL  ,NN1A|C2»C(.  »F,-?R 
1  FORMAT (5I=,hF1P. 5) 

31  C3=C2/C4 

C  N=  TOTAL  HUMBER  CF  FOINTS  READ  IN  a  LONG  BOOT. 

C  NSHA°F  is  A  PARAMETER  WHICH  DESCPICES  THE  30DY  SHAPE. 

C  NN1  =  NUMBER  OF  GRID  POINTS  COfPUTED  ALONG  first  0GIVE7NN2  ALONG  2 NO 

PORTION  OF  6CCY?  NN 3  ALONG  THIRD  PORTION  A  NO  NN  •»  ALONG  ATM  SEGMENT , 

C  MAXUMUM  CF  A  SEGMENTS  ALLOVABLE. 

N3  =  l  FOR  CONICAL  30ATTAIL,=?  FOR  OGIVAL  BDATT A IL .  IF  OGIVAL  BOAT T AIL 
IS  PRESENT  THFN  AT  LEAST  5  POINTS  MUST  BE  GIVFN  ALONG  BOATTAIL. 

C  C2  IS  A  FACTOR  WHICH  DETERMINES  STEP  SIZE  IN  X  DIRECTION. 

C 

C  POINTEC  BODY 

C 

C  NBlUNT=1 

C  C2=0 . 9  AND  C A=  20.  APE  NOMINAL  VALUES  FOR  THESE  PARAMETERS. 

C  NS H APE = 1 ?  NOSE  ONLY. 

C  NSHAPE=2?  NOSE  PLUS  AFTERBODY. 

C  NS HAPE  =  3?  NOSE  WITH  A  DISCONTINUITY  IN  IT.  THERE  MAY  OR  MAY  NOT  BE 
C  AN  AFTERBODY  PRESENT. 

C  NS H A°E  =  5?  NOSE  PlUS  AFTERBODY  PlUS  BCPTAIL. 

NS  HAPE  5?  NOSr  WITH  DISCONTINUITY  IN  IT  PLUS  AFTERBODY  PLUS  B0AT4IL. 
Nl=NUMOfR  OF  POINTS  ALONG  FIRST  00IV57N2  =  NUMBER  OF  POINTS  T  HP  OJ  G  H 
C  SFCOHO  OGIVE  INCLUDING  FIRST  OGIVE. 

IF  NSHARF  =  3  OR  5  ,  AT  LcAST  FIVF  POINTS  MUST  3c.  REAO  IN  ALONG 
C  EACH  OF  THE  OGIVES,  EVEN  IF  THE  OGIVE  IS  A  STRAIGHT  LINF. 

r ; 

C  HLUNTfd  BODY 

C 

C  N3LUNT  =  2 

C  C2- ,  05  AMO  C<*=  1.0  ARE  NOMINAL  VALUES  FCR  THESE  PARAMETERS. 

D  NFL  =1  FCR  SPHERICAL  CAP?  NB_=2  FOR  TRUNCATED  NOSE. 

C  WHEN  THF  BODY  IS  BLUNTED  NShAPf  MUST  BE  EITHER  S  OR  5. 

NSHAPf  -  3  , N‘ N 1 A s  1  ?  BLUNTED  NOSE  WITH  NC  DISCONTINUITIES  OThLF  THAN  THE 
C  INTERSECTION  OF  THF  CAB  with  OGI'/E. 

NSHAPFr  3,  N  N 1 A  =  ?  ?  3LUNTEO  NCSD  WITH  A  Ul SCONT I NU I TY  IN  THE  OGIVE  SO  THEKF 
C  ARE  ?  OGIVES  FRESENT. 

C  NSHAPE  =  F , NN1A  =  1?  R AMF  AS  A  )CVF  EXCrPT  FOATTCIL  PRESENT. 

*ISHAPC=5,NN1A  =  ??  SAME  AS  ArOVE  EXCEPT  BOAT  TAIL  P-'ESENT. 

IP  NN1A  =1  ,  T  H r  N  N 1  =  1  AIIg  N2.GF.P?  ]P  N)'1A  =  2,  THEN  N1.GE.5  AND  12.  Gc.9 
C  PR  =  RADII'S  CF  SPhECICAl  0A°  In  CALIBrFSIO0  TRUNCATED  PORTION). 

I  J  =  1 

IF( ICOUNT.GT.l)  GO  TO  32 

WRITf  (l  ,  *(.) 

'!•  FORMAT  (?(.X,*rODY  COO?  D  INATf  S  *  ,  //  ,  ?fc  X  ,  *  X  *  ,  1 1  X ,  *iR»  ,  /  ) 
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RAGE  2 


61 


70 


75 


80 


85 


90 


95 


100 


105 


110 


03  ?  1=1 ,N 
900(5,1)  X  ( I )  , 9  ( I ) 

W9IT*(*,33>  X  ( I )  ,9  ( I) 

33  F  03MA  T  (2DX,?F12.4) 

2  CONTINUE 

3  F09HA  T (2F15.1C) 

32  IF  (N3LUNT .  c0. 2)  CALL  BLUNT 
IF (NBLUNT .EQ.2)  GO  TO  5 
X  3  ( 1 >  =X ( 1 ) 

RR  (1)  =R(  1) 

IF (N1.NE.2)  GO  TO  4 
IF ( NN 1 A . EQ • 2 )  GO  TO  4 
510  00  508  1=2,5 

R9(I)=(9(T)-9(1))/(X(I)-X(1>) 

T  A  =RP ( I ) 

T  A8E=9ET  AM  A 

I F  (T  ABE  •  L  T  .  0 . 94  >  GO  TC  509 

508  CONTINUE 

509  XB(2)=X(I) 
f?B  (2)  =R  ( I ) 

9P<1> :PO (I) 

NN1  =  2 

K  =  1 

9BP(1)=9P<1> 

9BP(2)=9BP(1) 

GO  TO  5 

4  00  6  J=1 , 5 
L  =  1 

CALL  F0R5(X,R,X  (J)  ,RP(J)  ,N2,L) 

6  CONTINUE 
TA=9P  (1) 

TAdE=RFTA»TA 
RBP(l) =TA 

IP  (N1.CQ.2)  GO  TC  510 
TF(TA3E.LT..94>  GO  TO  503 
00  505  1=1,5 
TA9E=BETA*RP(I) 

IF  (TAB!-. LE. 0.94)  GO  TO  506 

505  CONTINUE 

506  X9  (2) =X ( I ) 

98 (2) =P( I ) 

9BP!2)=9P(I) 

RP(1J»PP(I1 

99(11=0. 

99P(1)=RP(1) 

XB(1)*X3(2)-PB(2)/9P(1) 

XB (3) =XB ( 2) *0 . 0 1 

J  J  =  3 
JK=? 

GO  TO  C0T 

5  03  CALL  FC5  <X(1),X(1),X(?),XC3),X<4),X(5),9‘M1>,RP<2>,9P<3),9F<4), 
1 9  P (5 ) ,B°F) 

PHOB=A!  S  (  «  1  .  *FP(  1)  »*?)  ♦»1.E/P.PP> 

XB(2)=0.0  2f  •‘•’HCB/9ETA**1.5  ♦XB(1) 

JJ-? 
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1  15 


1?0 


1?5 


no 


1  ?5 


190 


lL<i 


isa 


is=; 


1*0 


I'-'- 


a-  •"  "•'-•V*  r,f-  ■>;<ut  OPT='i  DR/i?/?/  16.20.37. 

J<=1 
JU  =  M 

ip (‘Mt.r..?) 

5  i. 7  J  -  l 
TO  T 

'All  ,V  j(K)  ,*8<K>  ,JU,3) 

CALI  FC°5  (X,  p,  y  3  (<)  ,PSD(K)  ,  J'J  ,  J) 

)£Tl=RfT5 

T  F  (BC?l , ST .1 . )  B?T!  =  l. 

X3(K*1  )  =  X“(<)  *P-T1  *  (R9U)  -i1)  (,!».))  *C2 
IF  (X^MK  +  D.GF.XIND)  GO  TO  3 
7  CONTINUE 
3  XT  <<*1 1  =  X (Nl  ) 

R° (K+l )=9  (Nl ) 

NN1=K*-1 
NN2=NN1  *13 

call.  fnP5(X,v,xa(K;  +  l)  ,P.JF(K  +  1),JU,  J) 

=  GO  TO(°fl'», 11  ,12,11)  , N3H& FF 

9  NN  =  N»Ji 

1 NL-X  0  INN) 

u=c. 

aia=o. 

GO  TO  90 

10  X9<K+2)=X (Nl) 

=*0(<+2)=R  (Nl) 

R0P(K+2) =0. 

BiTlzorTfl 

IF(9CTi.GT.l.P>  3£Tl= 1 , 

X<3(K  +  3)  =  C7  *3ETl*Rt)  (K  +  2)  ♦XO(K*2) 

R9(K*T)=RP(K*2) 

R0P(X*3> =0. 

H  <=<*! 

I J=T  J  *1 
a  - 1  j 
r 5r4»;T 
TTTi  =  ^ra 

I F  ( [lc  T  t .  G  T  .  1.0)  9c  T  2  =  1 . 

<ij  (K*')  =  Ce»°eTl*99  (K«?)  *X9  K-2) 

''?(K*’I)=:>P(<.;  1 
R5P(K*X)s0. 

(FIX.  "(•'*?). IT  ,  (  !  M  )  GO  TC  2  4 
X  P  ( K  *  X  )  =  X  ( M 
3f>  <<♦  <)»P.  (N) 

NN  =  K  + 7 
NHc’=X  *  ? 

ANL^X-1  (HM1) 
iL  =  C. 

a l  a=x 9 1’ri  >-xp  (mni  ) 

•»JT=HH2*  !  n 
GO  to  ci 

i ;  »"(<*: >  =  x (»  i ) 

)  =  9 (n; ) 

|>  M( 

"f-LL  r  P5(>  ,>-,xp(<,2)  ,->•■•<•  (•. .  V)  ,f 
F  ‘  ;  -  r  T  i 
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I F  CPrT 1 . 5T . 1 . 0 )  ?ET 1= 1 . 
XB(K+3)=X6(K+2)+C3  *  BET 1 *RB (K+2 ) 

CALL  INTERP  (X,P,XB(K  +  3) ,POJ<*3) ,K2, 3> 
CALL  FD°5 ( X  ,  R  ,  XB  (K  +  3) , R8F < K+3 > , N2 , J ) 
15  <»<+l 

IJ=IJ+1 
A  =  I J 
C  5  =  A*  C3 
BET 1  =  BFT  A 

XF(BETl.GT.l.O)  EETl»l. 

X9 (K  +  3>  =  XB  <K  +  2>+C5*8ETl  +  R8<K  +  2) 

IF1XB  CK+3) .GE.XIN2) )  XB ( K +  3 ) =X <N2 ) 
CALL  INTEPP (X,R,XB(K+3>,RB(K+3) ,N2,3) 
CALL  FCP5(X,R,XB(K+3) , RBP ( K  +  3 ) , N2 , J ) 
IFfXB (K  +  3) .LT. (X (N2)-.0001) )  GO  TO  15 
I  J  =  1 

X8 (K  +  3  >  =  X (N2  > 

RB(K+3)=R(N2) 

NN2=K ♦ 3 

ANL=XB(NN2)+RR 
IF (NFL.EQ.2)  ANL  =X  B ( NN2 • 

9L  =  0. 

ALA=0. 

CALL  FDP5(X,R,XB(K+3> ,*BP ( K  +  3 ) , N2 , J ) 

IF <X9 <K  +  3 > .LT. (X  IN)-.  000  1)  >  CO  TO  30 

NN=K+3 

NN3=K+3 

CO  TO  99 

30  RBP(K+4)=o. 

XB(K+4)=XB(K+3) 

RB <K+4)=RB (K+3> 

0ET1  =  BET  A 

IF(BETl.GT.l.O)  BET 1=1 . 

XB(K+5)=C3  *BETl+RB(K+4)+X81K+4) 

XB(K+5)=C3/lCO.*UETl*RB(K+4)+XB(K+4) 

RB(K+5)=RB<K+4> 

RRPIK+6) =0. 

16  K  =  K  + 1 

I J=IJ+1 
A  =  I J 
C5=A*C3 
BET 1  =  9FT  A 

IF (BETl. GT.  1.0)  BET 1=1. 
X9(K+5)=C5*BFTl*R8(K+4)+XR (K+4) 

RB (K  +  5 ) sR  B ( K+4 ) 

RBPIK+5) sO. 

IF(XB(K+5).LT.X(N2+1) )  GO  TO  16 
XB(K+5)=X (N2+ 1 ) 

RP(K+5)=R (N2+1 ) 

NN3=K  +  F 

A  L  A=X  B (NN3) -XB (NN2) 

IF(NSHAB£,E0.5)  GO  TO  13 

NN=K+5 

GO  TO  99 

12  X3(K  +  ?)  =  X  (N1 ) 
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225 


230 


235 


21.0 


21.5 


250 


255 


260 


265 


27r 


7'5 


1  7 


51’ 


20 


21 


9P(K+2)=:?  (N1  ) 

R8P  K*2)  =0. 

ANL=X"(NN1> 

'>FTl  =  8rTl 

IF(BETl.GT.l.O)  SET  1=1. 

Xd  (K+3)=C3  »PET1»R3(K»2)«XB(K*2) 

XB(K*3)=07/100.*PET1*RB(K*?)»XC(K*2) 

RB(K+3)=9B(K*2) 

9 BP  (K  *3 )  =0. 

K  =  K*1 
IJ=IJ*1 
A  =  IJ 
C  5  =  A*  C3 
BET1=BETA 

IFIBET1.GT.1. 0>  PET1=1. 

XB (K+3)=XB ( K  +  2 ) *C5*  BET 1*  RB (K+2J 
9fl(K+3)=9B(K*2) 

RBP(K*3>  =0. 

IF(XB(K+3) .LT.X(N1*1) )  GO  TO  17 
X9(K+3)=X (Nl»l) 

9B(K+3)=9(N1+1) 

NN2=K*3 

ALA=X9(N‘J2)-XP(NN1) 

IJ  =  1 

XB(K+4)=XB(K+3> 

RB ( K*4 )  =  PB  IK +3) 

R  BP  <K  +  4) =  RBP ( K ♦ 3  > 

BETl  =  eET  A 

IF(BET1.GT. 1.0)  eETl=l. 

XB(K*5)=XB(K  +  A)tC3  *  BET  1*9B  ( K+4 ) 

IF1N3.E0.2)  GO  TO  20 
XFlICOUNT.GT.l)  GO  TO  517 

SLOPE= (9  ( N) -PB ( K  +  4  >  ) / (X(N)-XB  <K  +  4> )  +1./57.293 
IF(  5LOPE.LT .-0.0372)  SLOPE=-. 1S72 
R(N)=SLOPF+ (X(N)-XB(«+4> ) +RB  (K  +  4> 

93(<*5)  =  9B(K+4)*- SLOPE*  (XB  ( K+ 5 > -XB ( NN2 > > 

99P (K+4) =SLOPE 
9.3P(K+E)  =SLOPE 
GO  TO  21 

CALL  I  NT E 95 (XP  (K  +  5) ,X (N-4) , X(N-3> , X (N-2)  ,X(N-1)  , X(N)  ,R(N-4>  , 

19(N-3)  ,9  (N-2)  ,9(N-1)  ,R(N>  tP.B(K  +  5>  > 

CALL  F05 (XP  OC+5 > ,X (N-A) , X (N-3) ,X (N-2) , x(H-l) ,X(N> , 9(N-A)  ,9 (N-3) , 
19IN-2) ,9(N-1> ,R(N> ,RBP(K+5)) 

IF(XR(<+  5)  .LT.X(Ki)  )  GO  TO  13 
X  3 (K  +  5 )  =  X  (N) 

R0  (K  +  5  )  =9  ( N  ) 

NN3=K+5 

3L=XB (NNJ)-XP (NN?) 

MN=K+5 
GO  TO  GO 
<:<*! 

I  J -I J  *  1 
A  =  I  J 
C5=A*C3 


0FTl=P '  TA 


19 
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325 
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22 


23 


13 


516 


24 


25 


IF  (BETl.GT.  1 .  (l)  SET  1=1. 

X3 (K*5)=XB(K«4) ♦  C5  *BE  T 1*R9  (<  *L ) 

I F  (N3  .  EO  .  2)  GO  TO  22 

R9 (K*5)=RB(NN2) *SLCPG*(XB  ) -X °  ( NN2 ) ) 

R3P (K  *r  >  =  SL  CPE 
GO  TO  23 

.3^I:LI^NIeR5(Xe,K  +  5,  »X  ‘N-*.)  ,X(N-3),X(N-2),X(N-l),X(N),R(N-4>  , 

1R (N-3) ,R (N-2) ,R(N-1> ,R(N) ,PB(K+5) ) 

““  ™5,XB'K*5,’X  (N-‘*>*x(N-3>.7<N-2)  ,X(N-1),X(N)  ,R(N-4)  ,R(N-3)  , 
1 R  ( N-2  )  ,R(N-1),R(K)  , R3  P ( K  »  5  >  >  ’ 

IF <XB (K*5) .IT. X (N) )  GO  TO  15 
XB(K+5) =X (N) 

R8(K*5)=R(N) 

NN3=K*5 

NN=K*5 

BL=XB (NN3) -XBINN2) 

GO  TO  99 
X8  (K*6)=XB(K+5) 

R8(<*6)=RP(K+5> 

IJ=1 

BET  1=  BFT  A 

IF(BETl.GT.l.O)  BET 1=1. 

XB(K  +  7)=XB(K'F6)+C3/2.*BET1*RB  (K*6) 

IF(N3.EQ.2>  GO  TO  24 
IF  (ICOUNT.GT.l )  GO  TO  516 

SUOPE= (R(N)-R3(K*6> >/ (X(N)-XB (K  +  6) )  ♦1./57.293 
IF (SLOPE. IT. -0.0 87 2)  SLOPE=-.0872 
R  CM) =SLOPE* (X(N)-XB(K+6) ) +R9 ( K*6 ) 

RB (K*7)=RB(K*6) ♦SLOPE •( XB ( K ♦ 7 ) -XB ( NN3 > > 

RBP ( K+6) =SLOPE 
RBP(K*7) sSLOPE 
GO  TO  25 

CALL  I NTER5 (XB (K*7) , X (N-4) ,X (N-3)  , X  (N-2 ) , X ( N- 1 )  , X ( N) , R (N-4 )  . 

1R (N-3)  iR (N-2) ,R(N-1) ,R(N) ,RB  (K‘7) ) 

CALL  F05 (XB (K*7)  ,X  (N-4) , X  (N-3) ,X(N-2)  ,X(N-1) ,X (N) ,R(N-4) ,R (N-3) , 
lR(N-2) jR(N-l) |R(N) *  RBP (K  *7) ) 

IF(XR ( K ♦ 7 ) .LT.X(N) )  GO  TO  19 
XB(K*7)=X (N) 

RB ( K+  7 )  =  R ( N ) 

NN4=K*7 
NN=K*  T 


BL-XB(NN4)-XB (NN3) 

GO  TO  99 
19  X=K+1 

IJ=IJ*1 
A  =  I J 
C5=A»C3 
BETl=  BET  A 

IF(BETl.GT.l.O)  BET1=1 , 

X3(<*7)  =  XB(K*6) *C5  *8GT 1 *RB  (K*6> 

IF (N3.EQ.2)  GO  TO  26 

RB(K*7>=RB(NN2)»SL0PE*(XB(K*7)-XB(NN3)> 

RBP(K«*> =SLOPE 
GO  TO  27 

26  CALL  INTERS IXP(K +7) ,X (N-4) ,X (N-3)  , X  (N-2 ) , X (N-l ) , X ( N> , R(N-4> , 
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SU940lJTINc  HIPPIC 

Cn“M?N/6£0“/'l-  (6)  ,X(30)  ,->(33)  ,r2,N,NSHA’E,Nl,N2,  X9  (22  5 )  ,  *  9  ( 225) 
COM^ON/GECl/  pjp  (225)  ,0>ETfl 

C04MON/GE02/NN1 ,NN2,NN3,NNt  ,NFL,  WOL  UNI ,  NNf  NNI , IHRINT.NN1A 
5  C0MMON/GFO3/0CtfS,AL ,XM,Y“, X  INI, Y INI, NNI ft 

CCuMON/GE04/K,f  ,=R,R-<FF 

C0«M0N/9IS2/  SUM1,SUM?,SUM3,SUM4,SUM5»3UM6 
COMMON/OAT1/  T(100),AK<liiD),aE(100>,C(225>,Cl(225>,C3 
COMMON/DISC/  I ,JK,AI2,SUM, JH,PI 
10  C0MM0N/DIS1/J1 , J3 

COMMON/WAvE/CADl  ,CNBL  » CMR L ) Cfl W » CN W , CMK 
C 0 MMCN/C  P V/  CPV  (225,20) ,JA,J9 

DIMENSION  PSI (’25) ,PMI (225) ,ZE0X (225) , Z E  OR  ( 225 ) , ZE 0 ( 22 5) , 

1ZE0XX (22  5), ZEOXR  (2  25)  ,ZE0RR( 225) , PSIX (225) ,P SIR ( 225) , ZtOP ( 22 5) , 
15  2ZE0PX (22  5) ,ZE OPR (2?F) , PHI  X (225), PHIR (225) ,0(225)  ,ZE1  (225)  , 

1ZE1X(225) ,ZE1R(225> 

DIMENSION  THET  (20)  •  I rt t T 1  (20 


DATA (T (I ) 

,1  = 

1 ,99)/. 01 

,.02 

,.03 

» ■ 

04, .05 

,.06, .07 

9  •  0  8  9  •  09 

,.10 

,.11,  .12 

1 

.13, .14,  . 

15, 

.16, 

.17,. 

18,. 

19,. 

2C 

1  •  2 1 ,  • 

22. 

.23,. 

24, 

.25,. 

26,  . 

27, 

.28, 

2 

.29, .30, . 

31, 

.32, 

.33,. 

34,. 

35,. 

36 

,.37,. 

38, 

.39,. 

40, 

.41,. 

42,. 

43, 

.44, 

3 

.  45, . 46, . 

47, 

.48, 

.49,. 

50,  . 

51,. 

52 

»  .53,. 

64, 

.55,. 

56, 

.57, . 

58,  . 

59, 

.60, 

4 

.61, .62, . 

63, 

.64, 

.65,. 

66,  . 

67,  . 

6e 

,  •  69,  . 

70, 

.71,. 

72, 

.73,. 

74,. 

75, 

.76, 

5 

.  77, .  78,. 

’9, 

.  80, 

.81  ,. 

82,. 

8  3 »  • 

84 

.  .85,. 

86, 

.37, . 

88, 

.89, . 

90,. 

91, 

.92, 

6 

.93,.  94,. 

95, 

.96, 

.97,. 

98,. 

99/ 

?5  DATA( AK( I) , 1=1,99) /3. 3590 2,3.02571,2. 63592,2.70215,2.6010 7, 

12.51987,2.45234,2. 39475,2.34473,2.30064,2.26132,  2.22592,  2.19380, 
22.16aa5,2.137a8,2. 1125’, 2 . 08946 , 2. 0679a , 2. 047 82 , 2. 02896,2. 01123, 
31. 9945 1, 1. 978 71, 1.963 76, 1.9-.957, 1.93608,1. 92324, 1.91099, 1.89929, 
41.88811,1.87740,1.86713,1.85727,1.84780,1.83870,1.82993,1.82148, 
30  51.81331,1  .8  0547,1.797  87, 1 .79053,  1 . 7 8343 , 1 . T7 655 , 1 .  769d9 ,  1.  76344, 

61.75718,1.75111,1.74521,1.73948,1.73392,1.72851,1.72324,1.71812, 
71 . 7131  3,1 .70  8  27,1. 70354,1. 69892,1 .69442 ,1.690  0  3,1. 68575,  1.68157, 
81 .67748, 1.67350,1.66960, 1.66579, 1.66206,1.65842, 1.65485, 1. 65137, 
91 .6479F,  1.64461,1. 64133,1.63813,  1 . 63499 , 1 . 63 1 91 , 1 . b2889,  1. 62593, 
35  a  1 .623  0  3, 1.62019,1. 61789,  1.61465,  1 . 611 96 , 1 . 60 9 32 , 1 . 60672,  1 . 6 0 41 8 , 

91  . 60168, 1.5992 2, 1.596 80, 1.59443, 1.592 10, 1.589 81,1.  58 755, 1 . 58  53-* , 
Cl. 58316, 1.58101, 1.57890, 1.57683, 1.5 74 79, 1.5 727 8/ 

D aTa (AE ( I), 1  =  1 , 99 >/l. a  28 36, 1.049 70, 1.0 08 35, 1.0 85 26, 1.10085, 
11.11541,1.12909,1. 14204, 1.15433,1.16606,1.17727, 1.18802, 1.19835, 
<•0  21.20828, 1.21 7»6, 1.227 11, 1.23604, 1.24469,1.25307, 1. 26119 , 1 . 2690 7, 

31.27672, 1.28416, 1. 29139,1.29843,1 .30523,1. 31196, 1. 31847, 1. 32482, 
41.33102,  1.33707, 1.3429“, 1.34875,1.35439,1.35991, 1. 36531, 

“I. 37059, 1.37575, 1.18082, 1.38577, 1.39063, 1.3 95 39, 1.40005, 1.40 46 3, 
61.40911,1.41351,1.41783,1.42207, 1.42623,1.43032,1.43433,1.43827, 
45  71 .44214, 1.44594,1.44968, 1 . 453 36 , 1 . 45697 , 1 . 460 53 , 1 . 46402 , 1 . 46746, 

81 .47085,1.47417,1. 47745,1.48068, 1 . 4 838 5 , 1 . 4 3 698 , 1 . 490 06, 1 . 4930 9, 
91. 496 07, 1.4990?, 1.501 92,1.5047 7, 1.50 759,1. 51036, 1.51 31 0,1. 51 579, 
ai. 51 84“, 1.521 07,1. “2366, 1.52621, 1.52872, 1.53121, 1. 53365,1.63607, 
91.53846 ,1.54081,1, c4 31 3 , 1 . 5 46 42 , 1 . 5476 9 , 1 . 54992 , 1.55213,1.55430, 
“>0  Cl  .5564  6, 1 .55  850,1.560  69,1  .“62  76, 1 . 56-.E  J  ,  1 . 56632 ,  1 .568 “2/ 

Pis’. 1516927 
T  (100 )  =  1  . 

aK(!00)=PI/?. 

AE  (100)  =  AK(100) 

55  C  THIS  SUBROUTINE  COMPUTF.S  1)6  SfcONU  ORDER  AXIAL  a  NO  FlUol 
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60 


65 


70 


75 


50 


55 


50 


95 


100 


105 


0°-r*  G^CSS  r l CM  Pi  RTu-' !AT  ICN  V c LOG  I” Y  COMPONENTS.  THESE 
COMPONENTS  ART  THtN  CTMpINtC  TO  YIELB  A  HYBRID  SOLUTION. 

I  <<=1 
I  K  =  1 

T  HFT  ( 1)=0. 

THFT1 (1 ) =0. 

00  47  IJ=2,19 

THET1  (IJ)-TriETl(IJ-l)  *10. 

THET< I J) =THET1 (IJ) 757.29553 
47  CONTINUE 
17  T  A=RP ( 1 ) 

IF(IPRINT.NE.l)  GO  TO  118 
WRITE ( 6 » 140  >  VOVS 

140  FORMAT (//, IX, “PRESSURE  COFF  FTCTrNTS  AT  M  r  *,F6.3,/7) 

WBITE(6,41> 

41  FORMAT(7X,1HX,10X,1MR,10X,5HOR/OX,7X,3HCP  ,/) 

118  TA2=TA»»? 

C  •  1)  =  T A2/SQRT  ll.-BETA**2*TA2> 

C  CONICAL  SOLUTION  ,  SUBSCRIPTS 
F 1 1=  ARSECH(BETA»TA) 

F22=  SORT ( 1 . -BET  A* *2* T A2) 

ZE0(1)=(F22-F11)»C(1) 

Z E OX  < 1 »  =-C  C 1 » *F11 
ZE0R(1)=C(1) *F22/TA 
ZEOXX (1) =-l./F22*C(l) 

ZE0XR<1)  =  1.7(F?2»TA)*CU> 

ZEORR ( 1) =-l./ (F22*TA2) *C(1) 

C  PARTICULAR  SOLUTION  AT  TIP 
1  =  1 

AN=1.2*VOVS**?/BETA**2 

PSIX(I)=VOVS»*2*(<ZEO  (I)  ♦  AN*  TA  * ZEO R C  I > > *ZF0X X ( I)  *ZE0X(I)*( 

lZEOXt I)«AN»TA  *ZE0XR(I) I -0.75»TA  *ZE OR  ( I ) * *2* ZEOXR ( I ) ) 
°SIR(I)=VCVS**2*((ZE0 (I) ♦  AN*  T  A  *Z EOfi < I ) > *ZE OXR < I > t ZE OX ( I ) »  (  < AN 

1*1.)*7F0R(I)*AN*TA  *ZE0RP(I))-0.25*Zl3R(I)**2* ( ZEOR ( I ) ♦ 3 . *T A 

2*  ZE9RR ( I ) ) ) 

C  C OMPL I MENT A° Y  SOLUTION  AT  TIP. 

C1(1)=TA*(TA* (1 . *Z  E  OX (1) ) -PSIR(l) ) 7F22 
A  °=C1 ( 1 ) 7C(1 > 

Zc0P( 1)  =  A3*  ZFO (1) 

ZEOPX ( 1) =A9*ZE0X  (1) 

ZEOPR (1) =AB*ZEPR(1 ) 

C  TOTAL  SOLUTION  AT  TI P=  PARTICULAR  PLUS  COMPLIMENTARY. 

PHIX( 1)=PSIX [ 1 >  +  7E0PX (1) 

PHIR(11  =  °SIR(1)+ZE0P3  Cl) 

OB=(l .*ZE0X (I) )**2*ZE0R(I>**2 

CP01=2./<1.4*VOVS**2)*m.*O.Z*VOVS**2»Il.-Q8)>**3.5  -  1.) 

QB=<1 . «°HIX  < I > )**Z*PHIR(I ) **Z 

CPV<l,l>=2./<1.4*VOVS»»Z>*( (1 .*0 .2* VOVS* *2* (l.-QB) ) **3.5-1 . ) 
CP0?=CFV(1,1) 

IF(IP-'INT.NF.l)  GO  TO  UP 

WRITE  <F,l?»  Xe(I)  ,RB(I)  ,R9MI)  ,CF02 

119  TFINN.FQ.?)  GO  TC  35 

C  F I  PST  CROr°  AXIAL  FLOW 

OC  7  I=2,NN 
SU"=0 . 
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115 


120 


125 


130 


135 


1<*0 


145 


150 


1F5 


160 


160 


36 


502 

62 

63 

66 

65 


6  ft 


I F ( I . lF, ?)  GO  TO  36 

J2=r-i 

00  ft  J=2 , J2 

XI=X9(J-l)-0fta»RB(J-l) 

TaO=BFTa»=n(i)/(yq (I) -XI) 

IF ( Tflu. 3E . 1 . 0 )  Tau=,309999 

,S^u,M:9CTfl*c,J,MX3,I,-XI,*,S0ST,i--TAu**2./TAu-Tau.aR5ECH,raun 

C  ONT I NU- 
JH=NN1+1 

IF(I.LE.JH)  GO  TC  36 
CALL  DISCI 
J=NN1 ♦ 1 

XI=XB(J-1)-BETA*RR(J-1) 

TAU=BETA*RB(I)/(X9(I) -XI) 

IFdAU.GE.l.)  Tau=0. 999999999 

SUM=SUP-BETA*C<J)*(X3(I)-XI)*(SQRT(1.-TAU**2>/TAU-TAU*ARSECHITAU)> 

J  H  =  NN  2  ♦  1 

IF(I.LF.JH)  GO  TO  36 

can  oisci 

J=NN2  *1 

XI=XB(J-i)-9ETA*R9<J-l) 

TAU=BETA*RB(I) /(X8(I> -XI) 

IF  (  Ta  U.GE  .  1.  )  iatl=0. 999999999 

SUH=SU0-9ET A*C (J)*(XQ(I)-XI)» (SQRT (1. -T  AU**2> /T  AU-T AU*ARSECH (TAJ ) ) 
JH=NN  3 ♦ 1 

IF(I.LE.JH)  GO  TC  36 

can  disci 

J®NN3* 1 

XI=XB ( J-l)-BETa»RB(J-l) 

TAU=9ETA*RR(I) /(XB(I)-XD 
IFdAU.GE.l.)  TAL=0 .999999999 

SUM=SUf'-BETA*C(J)*(XB(I)-XI)*(SQRT(l.-TAU**2J/TAU-TAU*ARSECH<TAU)) 
xi=xn (i-i)-BETa»oB(i-i) 

TAU=BETA*RB(I)/(X3(I)-XI) 

IF (TAU.SE. i . )  TAU=0. 999999999 
Taui=BETA*RB(I)/(XB(I) -X8(l) ♦PFTA*RB(l) ) 

IFITAU.GE.l.O)  TAU=. 999999 

TT=C( 1) *BETA*SGRT< l.-TAUl»*2) /  TAU1 

DEN=BETA*  (XE  (I) -XI) * (SORT ( 1. -T 80**2) /TAO-TAU*ARSECH(TAU)  ) 

c  C I ) =  (P3P(I)-TT-SUM)/0EN 

I F (IKK-2)  62,63,64 

JL=NN1 

GO  TO  65 

JL -NN2 

GO  TO  65 

JL=NN3 

IF(I.LE.JL)  GC  TC  ’ 

D=AQS (F9P (JL+1 ) -C9P (JL) ) 

IF  (D.Gr.  0.0001)  GO  TO  6d 
~ ( JL*1>  =C(JL) 

IKK=IKK+ 1 
GO  TO  7 

C (JL* 1 >  =  (R8PCJL*1>  -RJP  (JL)  )/  (RD°(JL*1>  ♦f?cTA> 

IKK  =  I<“'*t 
CONTINUE 
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C  1=2  I S  2ND  POINT  ON  SURFACE 
DO  9  I=2,NN 
SUM1  =  0 . 

SUM2=0. 

179  SUM3=0. 

SUML=  0. 

SUM5=0. 

SUM6=0. 

C  J=1  IS  CONICAL  SCLN.  WHICH  WILL  Bfc  fl DQED  IN  BELOW. 

175  HQ  to  J=2,I 

XXI=XB(I)*BETA»RB<J-l>-xn(J-l) 
tau=peta»rb(I)/xxi 
IF (TAU.GE. 1 . )  TAU=0. 999999999 
F1=ARSECH(TAU> 

180  F2=SORT (l.-TAU**2l 

SUH1=SUH1-C (J) *XXI**2*( (1 .♦0.5*T4U**2) ‘Fl  -1.5»F2> 

SUM2=SUM2-2,»C( J)*XXI* (F1-F2) 

SUM3=SUM3  *BETA*C(J)*XXI»(F2/TAU-TAU*F1> 

S UHL = SUM A -2 . *C (J)»F1 

185  SUH5= SUMS »2.» BET A*C( J)*F2/TAU 

SUM6=SUM6-BET A**2*C(J)*(F2/TAU**2  ♦  Fl) 

19  CONTINUE 
JH=NN1*1 

IF(I.LE.NNl)  GO  TO  18 

199  CALL  0ISC2 
J=NN1 +1 

XXI=X8(I>  *BETA*R?( J-l) -X9(J-1 ) 

TAU=BETA*RP(I)/XXI 
IF  (TAU.GE.l.)  TAU=0. 999999999 
195  F1=ARSECH (TAU) 

F2=SQRT(1.-TAU**2) 

SUM1  =  SUH1*C  (J)  *XXI*»2*  (d.*0.5*TAU*»2>  »F1  -1.5»F2I 
SUM2=SUH2t2.»C  (J)»XXI» (F1-F2) 

SUM3=SUM3-8ETA*C(J)*XXI* ( F2/T AU-TAU*F1 > 

200  SUM4=SUHL+2.*C(J)*F1 
SUM5=SUM5-2.*PETA»C (J)»F2/TAU 
SUM6=SUM6*PETA*‘»2*C (J) * (F2/TAU**?  +  Fl) 

J  H=NN  2  *1 

IF(I.LE.NN2)  GC  TO  13 
205  CALL  01  SC 2 

J  =NN2* 1 

X  XI=XB (I ) +BETA*R9( J-l) -XB (J-l) 

TAU=BETA*RB(I ) /XXI 
IF(TAU.jF,i.)  TCU=0 .99999999° 

210  F1=ARSECH (TAU) 

F2=SGRT ( 1 .-TflU*«2) 

SUM1=SUH1 tC (J) »XXI**2* ( ( 1 . *9 .5*1 AU**2! *F1  -1.5*F2) 

SUM2  =  SUH2*2.»C  ( J )  *  X  XI  *  (F1-F2) 

SUM 3  s  SUH3-PETA*C(J)*XXI*(F2/TAU-TAU*F1) 

215  SUM4=  SUM  L  +  2.*C(J)*F1 

SUMF  =  SUH5-2.*('ETA»r.  ( J  )  *F?/ TAU 
SUM6=SUH6*PETA«*?»C (J)*(F2/TAU*«2  +  Fl) 

JH=NN3+1 

IF(T.LF.NNJ)  f.C  TO  18 
220  CALL  DISC? 
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229 


230 


235 


2  A  0 


265 


250 


265 


260 


265 


27  0 


J=NN3*1 

XXI=X3(I>*S£TA*Rn(j-l)-xB(J-l> 

TAU=BET2’9P(I)/XXI 
IF (TAU.GE.l . >  TAu=0.999999‘99 
F1=ARSFCH (TAU) 
r2=SQRT(1.-TAU**2) 

SUM1=SUM1*C(J)*XXI**2*((1.+0.5*TAU**2)*F1  -1.5*F2> 
SUM2=SUM2*2.*C(J)*XXI*  (F1-F2) 

SUM3=SUM3-BETA»C <J)*XXI*(F2/TAU-TAU*F1) 

SUM4=SUM4*2.*C( J)*F1 
SUM5=SUM5-2.*BETA»C(J)*F2/TAU 
SUM6=SUM6*BETA**2*C(J)*(F2/TAU**2  *  FI) 

18  XIMX1=XB(I)-XB<1)  +BETA*R9<1) 

TAU=3FTA*RB(I)/XIMX1 
IF (TAU.GE.l.)  TAu=0. 999999099 
F 1  =  ARSECH  CTAU) 

F2=SQRT(1.-TAU**2> 

ZEO (I ) =SUM1«XIMX1» ( F2-F1 )  *C ( 1 ) 

ZE0X(I)=SUM2-F1*C(1) 

ZE0R( I)=SUM3+BFTA*F2/TAU*C (1 > 

Zc0XX (I) =-l , / (XIMX1*F2)*C (1) *SUM4 
ZEOXR(I)=0ETA/(XIMX1*TAU*F2) »C(1>  *SUM5 
ZEORR (I) =-eETA**2/(XIMXl*TAU»*2*F2) *C (1) +SUM6 
QB=(1 . fZEOX (I) ) **2  ♦  ZE0R(I>**2 

CP01=2./(1.4*VOVS**2)*((1.*0.2*VOVS**2*(1.-QB) >**3.5  -  1.) 

CPVII ,1) =CP01 
PHIX( I >  =  ZF0  X ( I > 

PHIP{ I )=ZE0R ( I > 

IF (NSHAPE.NE.4)  GO  TO  503 
I F ( I . GT . NN2 >  GO  TO  505 
GO  TO  504 

503  IF ( I, GT. NN3 >  GO  TO  505 

504  CONTINUE 

C  SFCONO  OROER  AXIAL  SOLUTION. 

C  A.  “ARTICULAR  SOLUTION 

A N= 1.2* VO  VS* *2 /BETA** 2 

PSI(I)=VOVS»*2*(ZEOX(I)*(ZEO<I)+AN*R3(I)*ZEOR(I))-0.25*RB(I> 

1*  ZE0R ( I ) *  *  3 ) 

PSIX(I)=V0VS**2*( (ZFO  (I)  +AN*RB( I )*ZE0R(I) )*ZE0XX(I)  +  ZE0X(I)*( 
1ZEOX(I)*AN*RB<I)*ZEOXR(I))-0.75»RB(I>*ZEORII>**2*ZEOXRCI>) 
PSIR(I)=VOVS**2*((ZEO(I>*AN*RB(I)*ZEOR(I))*ZEOXR(I)+ZEOX(I>*((AN 
l*l.)*ZFOR(I)*AN*R3(I)»ZEORR(I))-n.25*ZEOP(I)*»2*(ZFOR(I)t3.*RP(I) 
2* ZEORR  ( I ) ) ) 

C  3.  COMPLIMENTARY  SOLUTION 

20  SUM=0. 

IF(I.E0.2>  GO  TO  37 
J3=I-1 

00  12  J-2.J3 

XI=Xd (J-l )-?ETA*?B (J-l) 

T  AU=3CTA*RB(I)/(X3(I) -XI) 

IF( TAU.GE .1 . )  TAu=0. 99999959a 

SUM=SUM*"ETA*C1(J)  *1X3(1) -XI  )*(SQRT(1.-TAU**2)/TA'J 

l-TAU»ARSECHITAll)  ) 

12  CONTINUE 
JH=NN1*1 
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i£-a 


SJoeClITIPE 


*->yg-»  i  "i 
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RAGE 


2?3 


2»5 


29Q 


295 


30  3 


305 


310 


315 


320 


325 


330 


TRlCi 

Ic(I.Lb.JH  >  GO  13  32 
'ALL  nT9C3 
J=NN1*1 

X  1  =  X*?  (  J-11-0E  TAORB(J-I) 

TAU=0FTA*RB(I)/(X0<I)-XI) 

I F { TAU.5E . 1 • )  TAU=0.999999c9O 
SUM-9'JM-nETA»Cl  <  J>  *  < X0  <  I )  - XI  )  *  (SORT  (l.-TAU**2>  /TAU 

1-TAU*ARSSCH(TAU> ) 

JH=NN?*1 

IFd.L5.JH  )  CO  TO  .3? 

CALL  OISC7 
J  =  NN2* 1 

XI  =  XB<J-1)-BETA«SB  <J-1> 

TAU  =  0ETA*RBd)  /  <XB  (I)  -XI) 

IF(TAU.GE.l.)  TAU=0. 999999990 

SUM=SUW-BETA»C1(J)  MX0U1-XI  > * ( SQR T ( 1 . -T AU* » 2 ) /T A J 

1-TAU*ARSECH(TAU> ) 

JH=NN341 

IFd.LE.JH  )  GO  TO  37 
CALL  OISC3 
J=NN3+ 1 

XI  =  XB( J-1)-BETA*R9  (J-l) 

TAU-BETA’RB ( I ) /(XBCI) -XI) 

IF1TAU.GE.1.)  TAU=0. 999999999 

SUM=SUH-BETA*C1 ( J) * (XB(I)-XI  ) * (SQR 1 ( 1 . -T AU * * 2 ) /TAU 

1-TAU»A®SECH(TAU) ) 

37  TAU=BETA*RB(I) / (X8 ( I ) -X B < I -1 ) +PET A* R9 ( I - 1 ) ) 

I F ( TAU . G  E . 1 . >  TAU=0. 999999999 

TAUl  =  BETA»RB(I )/(XB(I )-XB( 1) ♦ BE  T  A»RB ( 1 ) ) 

T  T=C1 ( 1 ) *  BE  T  ASSORT <1.-TAU1'»2>/  TAUl 

Cl ( I)  =  (R9P (I) * ( 1 . *  ?E0  X (I ) ) -PS  IR  1 1 ) -TT-SUM) /  (BET  A*(XB(I)- 
lXBd-l)4BETA*RBd-l>)*<SQRT(l,-TAU**2>/TAU-TAU*ARSECH(TAL'))  ) 

I F ( IK-2)  99,95, 99 
99  JL=NN1 

GO  TO  92 

95  JL=NN2 
GO  TO  °7 

96  JL=NN3 

97  Ic (I.lE.JL)  GO  TO  03 

BxAoojprjp  (jl  *i)  -RBP(JL)  ) 

IFtO.LE. 0.0001)  GO  TO  92 
")  1  =  PS  I (JL41 ) -PSI ( JL) 

C  3=-01 

Cl ( JL  *1)  =  (P0P ( JL  +  1>* ( 1 .4ZECX (JL  +  l) >-PSIO (JL*1) -TT-SUM *3. *0  3/ 

1 (fl.'RB (JL  +  l) ) ) /OFT  A 
I K= IK ♦ 1 
GO  TO  °3 

92  CKJL41)  =  (Pnp(  jl  +  1)  *  (I.+ZECX  ( JL41  )  )  -P5IRIJL41)  -TT-SUMI/BETA 

I<  =  IK4  l 

9  3  SUM1  =  C . 

SUH2=0 . 

Sl)H3  =  0  . 

CO  13  J=?,I 

XXT=XJ<I) -XB (J-l) 4BFTA»Rfi( J-l ) 

T  AU=BFT 1  *  an ( I ) /XXI 
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345 


350 
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IFITAU.GE.l.)  TAU=C. 999999^90 
F1=ARSECH(TAU) 

F  2=SQRT (1 .-Tflu**2) 

SUM1=  SUM1 -C 1 (J)*XXI**2*<<l.*0.5*’TAU**2>,Fl-1.5*F2) 
SUM2=SUM2-2.*C1 ( J) »XXI*(F1-F2) 
SUM3=SUM3tOETA*Cl<J>'‘  ( F2/ T A U- TAU*F 1 > *  XXI 
13  CONTINUE 
J  H  =  NN 1  +  1 

IF(I.LF.NNl)  60  TO  21 
CALL  0ISC4 
J  =  NN1  + 1 

XXI=XB(I)-XB( J-1)*8ETA*RB<J-1> 

T  AU=8 F  TA  *  R8  ( I )  /XXI 
IFITAU.GE.l.)  TAU=Q. 999999999 
F 1=ARSECH (TAU) 

F2=SQRT(1.-TAU*«2) 

SUM1=SUM1*C1 ( J)*XXI**2» ( (l.+0.5*TAU**2)*Fl-1.5*F2) 
SUM2=SUM2t2.»Cl <J)*XXI* (F1-F2) 
SUM3=SUM3-BETA''C1(J)*  < F2 /T AU- TAU*F1 ) *  XX I 
J  H  =  NN2 ♦ 1 

IF(I.LE.NN2)  GO  TO  21 
CALL  0ISC4 

J=NN?+1 

XXI=XB(I)-XB(J-l)tBETA«R8(J-l) 

TAU=BETA*RB(I)/XXI 
F 1  =  ARSEC  H ( T  AU ) 

F2=SQRT(1.-TAU**2) 

SUM1=  SUN 1 *C 1 (J)*XXI**2*(  (l.+0.5*TAU**2)*Fl-1.5»F2) 
SUM2=SUM2n2.*C1(J)*XXI* (F1-F2) 

SUM3  =  SUM3-RETA*C1(J)MF2/TAU-TAU*F1),XXI 

JH=NN3+1 

IF (I.LE.NN3)  GO  TO  21 
CALL  DISC4 
J=NN3+ 1 

XXI=XB(I> -XB(J-1)+RETA*RB(J-1) 

TAU=BETA*Re< I) /XXI 
IF(TAU.GE.l.)  TAU=0. 999999999 
F1=ARSECH<TAU) 

F2=SORT<l.-TAU**2) 

SUM3=SUM3-BETA*C1( J)' ( F2/ T AU - T AU» FI ) *  XX I 
SUM2=SUM2+2.«Cl ( J) * XX  I* ( Fl-F 2 ) 

SUM1  =  SUM  1 +C1 (J) »XXI**2M  <1.  +  0.5*TAU**2) *F1-1.5*F2> 
21  XIMX1  =  X3( I) -XB( 1)  +8  ETA  *RP ( 1 ) 

T  AU=BE  T  A  *  RB ( I! /XIMX1 

IF ( TAU. GE .1 . )  T AU=0. 999999999 

F 1 =ARSECH {TAU) 

F2=SQRT(1.-TAU**2) 

ZE0P(I)=SUM1+Cl (1)*XIMX1*(F2-F1) 

2E0PX (I) =SUM?-C1 (1) *F1 
ZEOPR (I) =SUM7  +  C1  (1)*BETA*F2/TAU 
C  TOTAL  2N9  ORDER  SOLA. : PARTICULAR  PLUS  COMPLIMENTARY. 
°Hi(i)=psi(i) *  ZE  OP ( I ) 

PHIX(I)=PSIX(I)+ZEPPX(I) 

PHIR(I)=PSIP(l)+ZEOD^(i) 

98=(1  .-tPHTX  (  T  )  )  **2  ■»  PHIR(I)*»2 
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3-V(I,l)-?./(!.i.*\/pvS**?)*((l,»9.?'‘V?VS**?*(l.-19))**3,5-i.) 

1D0?=C f  V  (1,1) 

ft?  IFUPRINT.fF.  1 )  GO  TO  9 

WR ITE  (6 , 4 2)  X9 (11,43(1)  ,°BHI),CPV<I,1) 

42  F0RMAT(1X,LF1C.5) 

9  CONTINUE 

35  IF(ABS(au.GT.  0.001)  GO  TO  116 
10  117  1=1, NN 
on  117  j=i,iq 
c°vii,j)=cpv(i,i) 

117  CONTINUE 
GO  TO  113 

C  FIRST  O^OEF  CROSS  FLOW 
116  IF(IPRINT.NE.l)  GO  TO  120 
WRITE(6>49) 

49  FORMAT (//,30X,16H  HYBRID  SOLUTION,//) 

WRITE (6,51) 

51  FORMflT(7X,lHX,  10X,1HR,3X,6M  THCT A , 8X , 2HCP, //) 

120  J  =  1 

9(1)=2./9=TA/(SQRT ( 1 . -BET A * * 2 *T A2 ) / ( 3E T A ♦ • 2* TA2 ) ♦  ARSECh < BE  I  A *T A)  ) 

1  =  1 

T  A  U=BFT  A  *  TA 
F1=ARSECH (TAU) 

F2=SQRT(1.-TAU*»2) 

ZC1  (1 >  =9 (l>/2.  MF2/TAU-TAU*F1> 

ZE1X(1)=B(1)*F2/TAU 

ZE1R(1)=-BFTA*U(1)/2.»(F2/TAU**2*F1) 

00  53  IJ= 1, 19 

U3=C0S(AL)» (1. ♦PHIXd ) )*S2N(AL)»COS(THET ( I J) ) • ZE IX (I > 

V9=C0S(AL)*PHID(I) *SIN(AL)*COS(THET (IJ) ) *(1.+ZE1R(I) ) 

W9=*S IN ( A  L) *SIN(THET( IJ) >  *(1. +ZE1 (I) /TA) 

'19=U9»*2»  VB*  *  2  *W0»»2 

CPV(l,IJ)=2,/(1.4*VOtfS**2)*((l.*0.2*VCi)S**2*(l.-QrJ))*»3.5-l.) 

IF(IPRINT.NE.l)  GO  TO  53 

WRITE (6, 42)  XB(I),RB(I),THET1(IJ) ,CPV(1,IJ) 

53  CONTINUE 

I F (NN . NE. 2)  GO  TO  23 
00  131  I J=1 , 1 9 
C°V(2,IJ)=CPV(1,IJ) 

131  CONTINUE 
GO  TO  105 
23  J5=NN 

00  2?  1  =  2, J5 
SUM=0. 

J6=I-1 

no  14  J=1,J6 
IF(J.GT.l)  GO  TO  110 

TAU=RETA*RU(I) /(X3(I> -XB ( 1 > +BETA*RB t 1 ) ) 

IF  (TAU.GE.l. )  TAU=0.  )99999°9Q 
GO  TO  111 

110  TAU  =  9ETA*RB(I)/(X(l(I)-XB!J-l)*t>ETA,R3(J-l)  ) 

I F ( TAU • 3F • 1 .)  TAU=", 999999999 

111  rl=ARScGH(TAU> 

F  2=SQRT ( 1 .-TAU**2) 

IF(J.E".l)  GO  TO  1  "7 
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n=ABS (xi tj) -v b ( j-i ) ) 

TFID.LT. 0.000001)  GO  TO  14 
10’  SUM=SUM-9(J)*(F2/TAU»*2+F1) 

14  CONTINUE 

TAU=3ETA*R9 (I) / ( X 3 (I) -XB ( I  —  1 ) ♦9£TA*PB(I-1)  ) 

IF (TAU.GP.l . )  TAL=0. 999999999 
0  =  A3S (X9 ( I) -X9 (1-1 ) ) 

IF (O.LT. 0. 000001 )  GO  TO  114 

8(I)= (2 . /BET  A+SU8) / ( SCRT < 1 . - T AU* *  2) /TAU*» 2 tARSEC H < T AU) > 

60  TO  115 

114  8  ( I )  =  0 • 

115  SUM1=0. 

SUM2=  0 . 

SUM3=  0 • 

OO  15  J=1,I 

IF(J.GT.l)  GO  TO  112 

TAU=9ETA*RB(I> / (X9 (I) -XBCl) ♦PETA*R9(1> ) 

IFCTAU.GE.l.)  TAU=0. 999999999 
XXI=XB(I)-XB(1)  *BETA»RB<1) 

GO  TO  113 

112  XXI=XBCI)-XBCJ-1)+BETA*RB(J-1) 

TAU=BETA*RB(I)/XXI 
IFCTAU.GE.l.)  TAU=0. 999999999 

113  F1=ARSECH(TAU) 

F2=SQRT <1.-TAU**2) 

24  SUH1=SUM1»B  (J)/2.*  CF2/TAU-TAU*F1)  *XXI 
SUM2=BCJ)  'F2/TAU4-SUM2 
SUM3=-PETA/2.*BC J) *  C F 2/T A U** 2+F 1 > *SUM3 
IF(I.EC.l)  GO  TO  46 

15  CONTINUE 

46  ZE1 ( I ) =SUH1 
ZE1X(I)=SUM2 
ZE1RCI)=SUP3 
C  HYBRID  THEORY 

00  48  IJ= 1 , 19 

UB=COS  C  AL  >  * (1. +PHIXCI) l+SINCAL) *COS C THET C I J) > *Z E 1 X C I > 

VB=COS  (AL)»PHIR(I)  *SI N CAL >*COS( THET  (IJ)  )*  C1.*ZE1R<I>  ) 
H3=-SINCAL) *SIN(THET ( I J )  ) * (1 . +ZE1 C I ) /P9C I)  ) 
QB=UB**2+V9**2*WE»*2 

CPYtl, IJ) =  2./ (1.4»VOVS**2>  * ( C1.*0.2*VOVS**2*C1.-QB> >**3.5-1.  ) 
IFCIPRINT.NE.l)  GO  TO  49 

WRITE  (6,  4  2)  XB  (I  I,  RBI  II  .THETKIJ)  .CP<f(I,IJ) 

48  CONTINUE 

IFCN.EC.2)  GO  TO  27 
22  CONTINUE 

106  IF (NBLUNT.E0.2)  CALL  NEWT 
CALL  WAVE 
27  CONTINUE 
RETURN 
END 
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S  UORHUT I NE  INTfKMTX,  TY,X,Y,N,J> 

DIMENSION  TX  ( 1  r 0 ) , TY ( 1 00 ) 

I=C 

1  1  =  1*1 

5  IF (TX l I > .LE.X)  GO  TO  1 

IF(I.LE.J)  I  =  J 
IF (I.GT. (N-2>>  I=N-2 

CALL  INTERM  X  ,TX  (I-?),TX<I-1),TX1I>,TX(I*1),TX(I+2),TY<I-2),TY<I-1 
1>,TY(I),TY(I*1!,TY(I*2>,YJ 
RETURN 
END 
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SUBROUTINE  INTFRf<X,xi,X2,X3,X4,XE>F1,F2,F3,F4,F>5,F> 

C  F  POINT  LAGRANGE  INTERPOLATION  SUBROUTINE 

C  X1.LE.X.LE.X5 

a  1  =  cx-xe  > *(X-X3) *(x-xo)*  c  x-x5 » 

5  A2=(X-Xl)MX-X3)*(X-Xi.>*(X-X5) 

A3=(X-X1)  *{X-X2)*(X-X4)‘(X-XS> 

ao=(x-xi) *<x-x2>*<x-x3>',ix-x5> 

A5=(x-xn*(x-x2>*(x-x3)*(x-x4) 

D1=(X1-X2)*(X1-X3)  *  (X 1 -X4 >  *  (X1-X5) 

13  02=(X2-X1)» (X2-X3)* (X2-X4)* (X2-XS) 

03= (X3-X1) * (X3-X2) • (X3-X4) »<  X3-X5) 

D4=(X4-X1)  MX4-X2)  •  (X4-X3) * (X4-X5) 

05=(X5-X1) • (X5-X2) » (X5-X3)* (X5-X4) 

ClsAl/01 

15  C?=A2/D? 

C3=A3/D3 

,  C4=A4/04 

C5=A5/D6 

F  =  Cl*Fl  +  C2»F2  4-C3*F3  +  C4*F4  +  C5*FF 
20  RETURN 

ENO 
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SUBROUTINE  Nf WRAF(C7, RN,C6,F,CC1 ) 

C  THIS  SO  3  POU  TINE  USFS  NEWTCf  RAPHFON  “f  I  HOD  TO  SOH/£  FOR  MEAN 

C  SKIN  cu I C  T I ON  COtFFIClENT. 

CF=0.0025 
J  =  0 

1  F=C7/S03T (CF) -ALCG10 <RN*CF>  *CF 

13FOCF  =  -  •  5  *CT/  (CF**1.5)  -  »43P2°/CF 
J  =  J*1 
CF1=  CF 
CF=CF-F/OFOCF 
IFICF.LE. 0.0031)  CF=0.0001 
OCF=CF-CFl 

IF  (anS(OCF)  -1  .E-05)  2,2.1. 

4  IF  ( J-5  0 )  1,1,2 
2  CONTINUE 
RETURN 
END 
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SUDKOUTINE  NfWT 

COMMON/ GEOM/RP (6) , XI 30) ,=(30) ,02 , N , NSh apt, ni ,N2 , X9(22 5) , *9 (225) 
COMMON/GEOl/  R9P (225) ,  9FT  A 

COMMON/GE02/N.N1  , NN2.NN3,  NN4  , NFL,M3LUNT,NN,NNI,IPRINT,NN1A 

COMMON/GE03/VOVS,AL,XM,YM, XI  NT , YI NT , NN I  A 

COMMON/GE04/K.F  ,SR,R-.EF 

COMMON/ KAVF/CABL.CN8L , CMOL , C A W , CN W , CM W 

COMMON/CPV/  CPV(22E,20> , JA,J9 

DIMENSION  PH ( 2  0 ) , PH  1(20) 

PLPI=(1.2*VOVS»»2>  *»3.5* <6./(7.»V0VS**2-l.) >**2.5 

IF(NFL.EQ.l)  GO  TO  2 

CP0=  (0.90F*PLPI-1. ) / (0.7»VOVS**2) 

IF(IPRINT.NE.l)  GO  TO  19 
WRITE (6,6)  CPC 

6  FORMAT (1X.40HPPFSSURE  COEFFICIENT  ON  TRUNCATE0  NOSE  =,F10.5> 

19  CA=CP0’ (R (1 ) /RREF) 

CN  =  0. 

CM  =  0. 

CABL=C A 
CN8L=0. 

CM9L=0. 

XCP=0. 

CL=-CA*S IN ( AL) 

CD=CA*COS (AL) 

IF (NFL.E0.2)  GO  TO  20 

NNI=2 

GO  TO  99 

2  CP0=<PL»I-l.>/(0.7*VOVS**2> 

CS=COS (AL) 

SS=SIN(AL> 

PH(1)=0. 

°HI (1 >  =0 . 

IF  (AL.GT .  0.  0001  )  GO  TO  9 
NM  =  1 

GO  TO  10 

9  DO  3  1=2,19 

PHI (I >=PHI(I-1) .10 . 

PH(I)=PHI(I)/57, 29583 
3  CONTINUE 
NM=19 

10  IF  (IPRINT  .NE.l )  GO  TO  18 
WRITE (6, 8 ) 

8  FORMAT (IX, 39HP»ESSUPE  COEFFICIENTS  ON  S°HcPICAl  NOSE) 

18  X1=0. 

DX=(RR+XM  )/6. 

DO  4  1=1,7 
X2=X1-RR 

R?=S0RT(RR'*2-X2*»?) 

DO  11  L  =  1 , NM 
A=(1.-X1  /PR)*»2 

CP  =C°0*(A*CS**2* (XI  /PR-1.) *SQRT(1.-A)»COS(PH(L))*SIN(2.*Al) 
1  * ( 1 A ) *CCS (PH(L))**2*SS**2) 

IF (IPRINT .NE.l)  GO  TO  11 
M  R I T  F ( L , 5 )  X2, F2,»HI (L) ,CP 

11  CONTINUE 
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X1  =  X1 trx 
4  CONTINUE 

D  =  CP-CPV  <2,1°) 

n?=o 

no  12  1  =  3,  NN1 
X2=XB(I> 

IF1X2.GE.XINT)  GO  TO  15 
X1=RR+X? 

IF(Xl.GE.RR)  X 1=RR 
IF(Xl.GE.RR)  X2=RR 
R2=SQRT (RR**2-X2*»2) 

00  13  L=  1 , N fi 
A=(1.-X1/RR>**2 

CP  =C°0*(A*CS**2+(X1  /RR-1.)*SQRT(1.-A)*C0S(PH(L))*SIN(2.*AL) 

1* (l.-A)*COS <PH(L>) *»2*S3**2) 

IF(IPRINT.NE.l)  GO  TO  13 
WRITE (6,5)  X2,R2,PHI <L) ,CP 

13  CONTINUE 
IFO.GT.O.)  GO  TC  14 
01=02 

02=CP-CPV (1,19) 

IFIDP.LE.o.)  GO  TO  12 
SL0PE=(02-D1)/<XE(I)-XB<I-11) 

XNV=  XB(I-1)-D1/SL0PE 
NNI  =  I 
GO  TO  15 

14  01=02 
D2=CP-CPV (1,19) 

IF(D2. GE. 0. I  GO  TC  12 

SLOPE =(02-01) /(XB(I)-XB(I-l)) 

XNV=XB(I-l)-ri/SLOPE 
NNI=I 
GO  TO  15 
12  CONTINUE 

15  IF  (I.GE.NN1)  XNV  =  X INT 
NNI  =  I 

IF (X2 .GE. XINT)  XNV=XI NT 
IF(I.GE.NNl)  NN I =1-1 
YNV=SQRT (RR**2-XNV**2) 

TH2=ATAN (-YNV/XNV) 

S«=SIN (TK2) 

CH=C0S(TH2) 

RA= (RP/RREF) **? 

C  A=CP  0/2. *RA* (CS**2*(1.-CH**4)+.5*SS  **2»SH**4) 

CN=CP0*RA*SIN(2.*AL)*SH**4/4 . 

CM=-CP0/2.*RA*SIM2.,'AL) ,(SH**4/4.+Sh**2*CH*»3/5.*2./15.»(CH**3 

l-l.)) 

CAJL=CA 
C  NBL=  CN 
CN9L=CM 

Cl  =  CN*CS-CA»S<: 

C.)=CA»CS*CN*SS 
X  CP=-CM/CN 
70  CONTINUE 

5  FORMAT (IX, 46 10 .5) 
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SURRGUTIN;  NOR“FC 

COMMON/G£OM/°b,6),x(  >0) ,R (20) ,02,N,NSMAPE,N1,N2,XB(225> ,R9(225> 

C  OMMON/G  E  01 /  REP  (225)  , SETH 

C0MM0N/GP02/NN1 ,KN2,NN3,NN4,NFL,N9LUNT,NN,NNT, IPRINT.NN1A 

CO“MCN/GEO3/V0VS,AL,XM,YM, XINT,YINT,NNIfl 

COMMON/GEC4/K,)  ,PR,P*>EF 

COMMON/ WAvE/CABL,CN8L , CMDL , CA W, CN h , CMW 

COMmon/UOL/  VOL , Cfl F, CNF, CMF.RN, O I A , XP , A P , VOLN 

COMMON /LENG/BL  ,  ANL, ALA 

DIMENSION  A1(10>,A2(10>,AM(10>,F1(10>,G1(10>,F2(10>,G2{10>,FA(10) 
1,004(10) ,006  (  10) ,0  07(  1  0) ,D0  8  (10) ,  009  (10)  ,010  (  10)  ,D12( 101 
2,011(10) , XCPLfl (10) 

OflTfl(fll (I >,1=1,10) /l. 75, 1.82, 1.9, 1.06,2. 05, 2. 6, 3.5, 3.65, 3.7,3.35/ 
DATA (42 (I), 1  =  1 ,101/ 1.9, 1.93, 1.69, 1.95, 1.97,2. 15, 2. 45, 2. 44,  2. 4, 2.2/ 
OATA(AM(I),Irl,10)/0.,.2,.t,.6,.68,.9,.94,.97,1.05,1.2/ 
3ATfl(Fl(I),I=lflO)/0.,.l,.2,.3,.45,.6,.75, .85,. 925,1./ 

DATA(G1( I), 1=1, 10)/ 3. 35, 3. 48, 3. 6,3.65, 3. 5,2. 6, 1.75, 1.46, 1.35,1.28/ 
D ATfl(F2( I) ,1  =  1,10) /O. , .25, .5,  .65, .82,1. ,1.2  5,  1.5,2. ,2.5/ 

QATA( G2( I) ,1=1,10) /3. 35,3., 2. 5f ,2.25, 1.93, 1.3,1. 02,. 95,. 85, .75/ 

DA T A ( FA ( I ) ,1=1 ,10) /O. ,.5,1. ,1.5,2.  ,2.5,3. ,4. ,6. ,8./ 
DATA(DC4(I),I=1,10)/0.,0.,G.,G.,0.,".,0.,0.,0.,0./ 

OATA(0 06 (I), 1=1, 101/0. ,.03, .043, .05, .05, .05, .05, .06, .05, .05/ 

OATfl( 007(1), 1=1, 10>/0. ,.03, .113, .133, .143, .148, .15, .15, .15,. 15/ 
DATfi(D09(I)  ,1  =  1,10) /0 . ,115, .16,. 18b, .207  ,.223, .235, .248, .25,  .252/ 
DATA(D10(I),I=1,10)/0.,.175,.23,.265,.293,.31,.325,.337,. 34, .34/ 
DATA (01 1(1)  ,1  =  1, 10  I/O.,. 097,. 138,.  16,. 176,. 136,.  19,. 195,. 19 7,.  197/ 
DATA(D12 (I) , 1  =  1,10) /O. ,.  097, .138, .16, .176, .136, .19, .195, .197  ,.197/ 
DATA(XCPLA(I>  ,  1  =  1,10) /. 5,. 4, .342, .31,. 29,. 272,. 26,. 248,. 245,  .245/ 

I F  (BL • GE. 0. 02 )  GO  TO  6 
CNAIO=0. 

GO  TO  2 

6  IF  (VOVS.GT.l.)  GO  TO  1 
F 1 1=SURT ( 1 . -VO VS**  2 ) 

CALL  INTFRP(F1,G1,F11,G11,10,3) 

CNALB=-G11* (1 ,-4.*RB(NN) »»2) 

GO  TO  2 

1  F12=SQRT (V0VS»*2-1.) 

CALL  INTERF(F2,G2,F12,G12,9,3) 

CNALB=-G12*(1.-4.*RB(NN)**2) 

2  T  HE=A  BS  (RBP  ( NM ) ) 

IF(NN1A.E0.2)  THE=ABS (R0P(NN2> ) 

I F  (NBL'JNT  .E  Q  •  2 )  THE  =  A  BS  (R9P  (NM2)  ) 

CALL  INTERP(AM,A1,\/0VS,A11,1D,3) 

CALL  INTERP(AM,A2,V0VS,A22,1C»3) 

CNALN=-flll»THE  *422 
IF(ALA.GT.O.Ol)  GO  TO  9 
CN ALA  =  0 . 

GO  TO  4 

9  CALL  INTERP(FA,D07,ALA,0071,10,3) 

CALL  INTFRP (FA, 003 , ALA, 0031 ,10,3) 

CALL  INTERP( FA, DIO, ALA, 0131,10,3) 

CALL  INTFRP(FA,D36,ALA,0061,10,3) 

IF(VCVS.GE. 0.8)  GO  TO  5 
9041=0. 

CALL  I NTER5 ( VCVS , • 4 , . 6 , • 7 , . 9 , 1 . , D 041 , 0 06 1, 007 1 , U 08 1 ,010 1 , CNA LA) 
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SUBROUTINE  NCRHFO  TRACC 


CDC  6609  FTN  V3.U-P3U8  OPT  =  ii  19/12/7,;  16.20.37. 


60 


65 


70 


75 


IF(VOVS.LT.0.6)  CNALA=5.*OCSl'<vnvS-9.4> 

I  F  ( VOVS  .  L  T  .  0 . 4  )  CNALA  =  0  • 

GO  TO  L 

5  CALL  I NTE RF  (F  A ,  012  »  AL  A  , 0 12 1 , 1 0 »  3) 

0111=0121 

CALL  INTER5(V0VS,.7,. A ,1. ,1.1, 1.2,0071,0081,0101,0111, 0121,CNALA> 
4  CNAL=CNALA*CNALN*CNALB 

XCPN=ANL- VCLN/ (3 ,14159*RREF**2> 

XCPB=XP(NN) -°L /2  • 

XCPB=BL* < 1.-3. 1 4159* (PREF**2-RREF*6L*A9S  (RBP(NN) )*BL**2/3.* 
1R8P(NN)**2))*ANL*ALA 
XCPA=ANL*ALA/2. 

CALL  INTERP(FA,XCPLA, AL A , XCP1 , 10 , 3 ) 

XCP=XCP1*ALA 

XCPA=ANL+XCP 

CMAL=- (CNALN*XCPN*CNALA»XCFA+CNALR*XCPB> 

CNW=CNAL* AL 
CMW=CMAL*AL 
RETURN 
END 
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SUBROUTINE  R9ANr>  TRACE 


COC  6f 00  FTN  V3.Q-P308  OPT  =  1  09/12/72  16.20.07. 


SUIRO'ITI  NE  R^ANO 
COMMON/ 3SMO/C A P.CNP, C“P,  Hi 
COMMON/GE03/VOVS,A|_,XM,Y.M,XINT,YINT,NM4 
DIMENSION  4M(15),3CAP(15) 

5  OATA(A-m,I=l,l?)/0.,.5,.7,.8,.9,.95,l.l,1.2,i.5,2.,2.5,3.0/ 

OATA(UCAP(l > ,1=1,12) /0. ,0. ,0. ,.001,. 0li5,.01,. 0092, . 0079, .0067, 
1  .  0055, .0052, .005/ 

FALL  INTERP ( A M, OCA P, VO VS, CAP  I ,12,3) 

C  AP=C  AP 1  *  H8 /0 .01 
10  CNP=0. 

CMP=  0 . 

RETURN 

END 
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SUBROUTINE  SI  "P 


trace 


cor  6600  FTN  V3.0-R309  CF-T  =  G  09/12/72  16.20.07 


PAGE 
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SUBROUTINE  SI“F 

CONMON/G? CM/PF (6) . X ( 30 ) , R ( 30 ) , C2 . N , NShAPE, N1 ,N2 , X9 (225 ) ,  PB  <225) 

COMHON/GEOl/  R9P(225) ,9ETA 

CO-MOM/CPV/  CPV(225,20> , JA,J8 

C04M0N/0IS2/SL*1  ,SUM2 ,SU*3 ,SUM4, SUM5, SUM6 

C0MM0N/GEC4/K,F,RR,RREF 

DIMENSION  FI (225) ,G(225) ,  Gl (225) 

DO  1  I=JA,JB 

al=  CPV  ( 1, 1) *2.* (CPV (I,3)*CPV (1 ,5) *CPV(1 ,7) ♦CPV(I,91 

i+CPV<I,ll>*CPV(I,13!+CPV(I,15)*CPV(I,17) )  ♦CPV(I,19> 

A2=  4.' (CPV (1,2) ♦CPV(I ,4)*CPV(I,6) ♦CPV(I,8)*CPV(I,10> 

1*CPV(I,12)*CPV(I  ,14) +CPV( 1,16) *CPV( 1,18) > 
F1(I)=G.15818*(A1*A2)*RB(I) 

Bl=  CPV (1,1) -CPV (I ,19)  *2 . *  (.9  3969*  (CPV ( 1 , 3 ) -CPV ( 1 , 17 ) )  +  . 766 04* ( 

1CPV(I,5)-CPV(I,15))*.5*(CPV(I,7)-CPV(I, 13))*. 17365* (CPV ( I ,9) - 

2CPV (I , 11 ) )> 

9  2=4. *(. 98481* (CFV(I,2)-CPV( 1 ,18) ) ♦. 8660 3* (CPV (I , 4 ) -CPV ( I , 16 ) ) ♦ 
1.64279* (CPV (1 ,6) -CPV (1, 14) > ♦  .  342 0 2* (C PV ( 1 , 8 > -CPV  ( 1 , 12 >  > ) 

G(  I)  =  0.0581 8* (Bl*B2)*RB(I) 

Gl  (I)=G(I)*XB(I) 

1  CONTINUE 

IF(JA.NS.JR)  60  TO  2 
SUM1=  0 • 

sum2=o. 

SUM3=  0 • 

GO  TO  99 
2  J9B=JB-1 

DO  3  I=JA,JB8 
7  H=(RB(I*l)-RB(I))/6. 

X12=(XP( 1*1) +XB (I) )/2. 

IF((JB-JA).LT.5)  GO  TO  4 
J  =  JA«  Z 

CALL  INTERPtXB, FI, X12,F12,J9,J) 

CALL  INTERP(XP,G,X12,G12,JB,J) 

CALL  I NT  ERF ( XB , G 1 , X 12 , Gl 1 2 , J B , J) 

GO  TO  f 

4  F12=(F1(I)*F1 (1*1) 1/2. 

G12=(G(I)+G(I+l))/2. 

G112=(Gl(I)+Gl(I*l))/2. 

IF ( JB . GT .2)  GO  TO  5 
F1(3)=F1(2) 

G  ( 3)  =  G  (2  ) 

G1(2)=2,/3.’G1(2) 

Gl ( 3 ) =G1 ( 2) 

G112=2./3.*G112 

5  SUM1=SUM1+H*(F1 (I) *U. *F12*F1 (1*1) ) 

H1=(XB(I*1) -XB  (I)) /6. 

SUM2=SUM2tMl* (G(I) ♦4.*G12*G(I*1) ) 

SUM3=SUM3+H1» (Gl (I ) *4 . *G 1 1 ?♦ G 1 ( I *1 ) ) 

3  CONTINUE 

SU*3=SUM3*SUM2*RR 
9°  RETURN 

END 


CD 


■C* 

in 


S'JRPOUT  I NE  SKINF  TRACE 


C^r  66  0  G  FTN  Vi.u-Bj  u8  OPTrO  04/12/7..  lb.20.J7. 


5 


10 


15 


?0 


?5 


30 


35 


LQ 


45 


50 


55 


4 

5 


6 

C 

C 


SUBPOUT I  NT  SK  INF 

COf'«ON/GEOM/PP(6),X(33)  ,  -  (30)  ,C2  ,  N,  NSHAPE ,  HI ,  N2 ,  XFU225)  ,  »  3  (2  26  ) 
00«“0N/GE01/  REP (225) . nrT a 

COM“0 M/GE02/NN1 ,  NN2 , NN3 > NN4 , NFL , NBL UNT , NN , NNI , IPRINT , NN1 A 
C0MM0N/GE03/V0VS,AL,XM, YM,XINT,YINT,NNIA 
COMMON/  GE  O1*  /  K,F,RR,-iR£F 

COMMON/ OI S2/  SUM1,SUM2,SUM3,SUH4, SUM5.SUM6 

COMMON /CPV/  CPV(225,20) , JA.JB 

COMMON/VOL/  VCL,CAF,CNF,CMF,RN,JIA,XP,AP,VOLN 

IF(NBLUNT.FO.l)  GO  TO  5 

IF  (NFL.EQ.2)  GO  TO  5 

SUM1=6.28318»PR*YINT 

SUM2=3. 14159' Y I  NT* *2* (PR-YINT/7.  ) 

THE=ATAN(-YINT/XINT) 

SUM7=RP**2*THE/2 .“PR* *2* SIN (THE) /2. 
ag=ACOS( (RR+XINT) /RR) 

S'JM4  =  SUM3*2./3.*FR**3*  (l.-SIN(AB>**3> 

GO  TO  6 
SU“1=0. 

SUM2- 0 ■ 

SUM3=0. 

SJM4=0. 

CF3=0. 

THIS  SUBROUTINE  CALCULATES  the  axial  force  coefficient  oue  to  skin 
FRICTION  ON  THE  GOOY  (CDF). 

PI=3. 14159 
AREF=PI*RREF*»? 

GAMA=1 .4 

T WOTI=l. +0.9* (GAMA-1. )*(VDVS**2/2. I 

c =g am  a-i . 

A  =  SQRT  <C*VOVS**2/(2.*TWOTI) ) 

S= (1. + .5*C*VCVS*»2) /TMOTI-1. 

0=SQRT (3**2*4.*a*«2) 

Cl=(2.*A**2-tJ>/0 

D2=8/n 

D 3=.2 42/ ( A’SORT  (THOTI ) ) 

04=ASIN(C1) 

C5=aSIN(02) 

C6=(l.+2.*.76)/2.*ALOG10(TWOTI) 

C7=03* ( 94+C5 ) 

RN3=PN*XB (NN) *OI A 

CAuL  NFWRAP ( C 7 , RN3 , C6 , H , CF 3 ) 

TF  (NSLUNT  ,EO.  1)  GO  TO  1 

< - NNI A 

K1=NN1 

K2=NN2 

K  3  -  NN  3 

K  4  =  NN4 
GO  TO  2 
K  =  1 

Kl=NMi 
K  2  =  NN  2 
K  3=NN  3 
K4=NN4 
JA  =  K 
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3U8RCUT INS  SKIMr 


1RAC;: 


car  bf-Qo  ftn  v3.o-P3»,e  opt=o  1*9/12/72  ie. 20.07 


PAGE 


jq=Ki 

IF  (J8.E.3.NN)  JP  =  NM-1 
CALL  TRAPF 

VOLN=  SU!12 

60  IF(NNl.cQ.NN)  GO  TC  99 

JA=K1*1 
jn=K2 

IF(JB.FQ.NN)  J  B  =  AN- 1 
CALL  TRAPE 

65  I F ( NBLUNT .EQ.2)  V0LN=SUM2 

IF(NNJA,£Q,2)  V0LN=SUM2 
IF (NN2.EQ.NN)  GC  TC  99 
J  A  =  K?  *  1 
J  9  =  K3 

70  IF(JB.FQ.NN)  J  8  =  K'N  - 1 

CALL  TRAPE 

■  IF (NN3. EQ.NN)  GC  TO  99 

JA=K3+1 
JB=NN-i 

75  CALL  TRAPE 

99  SB  =SUM1 
VOL=SUM2 
AP=SUM3 
XP=SUMA/SUM3 

HO  COFB=CF3*S8/AREF 

CAF=COF9 
CNF=0. 

CMF=C . 

RETURN 

85  EN9 


* 


B-47 


SU3°CUTINE  TRANS  T»BCE 


C1C  660  n  FTN  V3.0-P308  OPT  =  0  09/12/72  lb. 20. 07. 


PAGE  1 


SU9PCUTINE  TUANS 

COM"ON/',ECM/RF(6)  ,X(30),R(30>  ,C  2 ,  N ,  NSH  A  PF  ,  Nl ,  N2 ,  X  3  (225)  ,  R9  (225) 

C  OMMON/GE  01/  PUP  (225)  ,  BETA 

C0MMCN/GE02/NN1,NN2,NN3,NN<.,NFL,NBLUNT  , NN, NN I , IPRI NT , NN1 A 
5  OOMMON/GE03/VOVS,AL ,XH,Y», XI  NT, TINT ,  NNIA 

COMMCN/GEO<-/K,F,RR,RREF 
COMMON/CPV/  CPVC225.20) , JA,J0 
COMMON/ 3 IS2/  SUM  1, SUM 2, SUM 3, SUM A, SUMS, SU«6 
COMMON/LENG/BL.ANL  ,AL A 

10  COMMON/WAVE/CABL,CNBL,CM9L,CAW,CNW,CMW 

9IMENSI0N  AM (10)  ,CA15(10) , C A 2 ( 1 0 ) , C A 2 5 ( 10 )  , CA3 ( 1 0 )  , CA4 ( 1 0 ) 

9ATA( A* (I), 1=1, 8)/. 85, .9, .95,1. ,1.05, 1.1, 1.15, 1.2/ 

0ATA(CA15(I),I=1 ,8)/.01,.072,.13,.177,.215,.247,.277,.3/ 
DATA(CA2(I),I=1,8)/0.,.036,.073,.107,.14,.!69,.191,.205/ 

15  0ATA(CA25(I) ,  1  =  1  ,8 > /0 . , .  0 1 , .  04 , . 0 7, . 098, . 122 , .13 8  ,  .  143/ 

OATA(CA3(I),I=1,8)/0.,0.,.024,.048,.073,.092,.102,.097/ 

D A  T  A ( CA4 (I) ,1  =  1, “)/0. ,0. , • 01 , .032, .047, .055, .055, .04/ 

IF (RBP(NN) .LE.  0. )  GO  TO  87 
C  BO=0 . 

20  GO  TO  88 

87  tf OV=VOVS- .0498 

IFIVOV.LT.l.)  GO  TO  89 
ARfcF=3.14159*RREF**2 
IF (NSHAPE.NE.4)  GO  TO  1 
25  IF (NN1A.E0.2)  GC  TO  1 

J=NN2*1 
GO  TO  2 
1  J=NN3+1 

2  DO  10  L=J,NN 

30  X  X  =  XB ( L ) -XB(J) 

DELTA=ATAN(1./(2.*ANL) ) 

IF(RBP(J-3> .LT.RBP(l) )  OELTA=ATAN ( . 2/AHL) 

GAMA=1 . 4 
C 1=1. +GAMA 

35  C  0  =  SORT ( Cl > 

C3=VOV  “2 
C  4  =  1 . -C3 
C5=C4/(C1*C3) 

C6=3.*DELTA/(2,*C0) 

40  C7=25.»C1*V0V  **(2./3.) 

C8=.5*C4/(C1*C3) 

C  9  =  1 . ?5*C5**2 

C 10  =  2 • *C  5/ ( VCV  •* (2./ 3. ) ) *C6* • (2./3. ) 

C11=(C6/V0V  )»»<4./3.) 

45  CSQ=C7»(C8*(C9*C10*C11>**  (.5)  ) 

C=SQRT (CSQ) 

T=2.*ALA  *2. • XX 

CPI  =.4* (Y-C)/SQRT(Cl»VOV  * * ( 2 . / 3 . ) ) * ( . 04* ( Y-C ) * * 2/ ( C 1  * VOV  ** 

K2./3.) )  -C4/(Cl*C3)  )»».5 
50  IF (Y.GT.C)  CPl=0. 

0ELTA=-R9P(L) 

GAMA=1.4 
C 1 = 1 . ♦  G  A  M  A 
C0=SORT(Cl) 

55  0 3  =  VOV  '*2 


09 

I 

i- 

00 


SUORO'JT  IN:  TRANS  1RACF 


COC  6600  FTN  V3.0-P3J8  OPT=0  09/12/7,: 


60 


65 


70 


T5 


80 


85 


50 


95 


100 


C  4  =  1 . -C3 
C5=Ca/ (Cl*C3) 

C6=3.*DELTA/ (2.*C1) 

C  7=?5 . *C 1  *  VO  V  **  (2./3.) 

C3=.5*C4/(C1*C3) 

C9=l. 25*C5* • 2 

ClQ=2.*C5/<V0V  •*<2./3.>),C6*M2./3.) 

CUMC6/V0V  )  **  (4./3.  ) 

CSQ=C7* (C8*(C9+C10*C11)»* (.5) ) 

O-SQRT (CSO) 
r=xx*2. 

CPVIL,  1>=.4*(Y-C)/SURT<C1*V0V  **<2./3.))*<.04»(Y-C)**2/(Cl»VOV  ** 
l(2./3.) >  -  C4/  ( C 1*C3 )  )  **.5-DtLTA**?*CPl 
IF(Y.GT.C)  CPV(L,1)=CP1 
IF(IPRTNT.NE.l)  GO  To  15 
WRITE (6. 13)  Xe(L),CPV(L,  1),CP1 
13  FORHAT(1X,3F10.5> 

15  DO  10  *1=1,19 
CPV(L,K1)=CPV<L,1> 

10  CONTINUE 
JA=J 
JB=NN 
SUM1=0 . 

SUM2=0. 

SUM3=0 . 

CALL  SINP 

CB0=2.*SUM1/AREF 

V0=V0V-1. 

IF(VO.LE.0.06>  Cei=C90 
V1=V0VS 

89  CONTINUE 
IF(VOVS.GT.0.95)  GO  TO  90 
CBO=0. 

GO  TO  88 

90  CBO=C9l* (VOVS-.95)/  (V1-.95) 

88  CONTINUE 

CALL  INTERP(AM,CA15,YOVS,AO,8,3) 

CALL  INTERP (AM.CA2  , VOVS , A  1 , 8 , 3) 

CALL  INTERP(AM,CA25,V0VS,A2,8,3) 

CALL  INTER°<AM,CA3  , VO  VS  ,  A  3 , 8 , 3 ) 

CALL  I  NT  ERP ( AM ,CA4  , VO  VS , A  4 , 8 , 3 ) 

IF(ANL.LE.4.)  GC  TO  16 
CAN=A4*(l.-.2*(ANL-4.) ) 

GO  TO  17 

16  CALL  INTER5 (ANL,1.5,2. ,2.5,3. , 4 . , A  0 , A l , A2 , A3 , A4 , CAN) 

17  CAM=CAN*CBO 
99  RETURN 

ENO 


16.20.37. 
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SUBROUTINE  TRAPE  TRACE 


COC  650  0  FT  N  V3.0-P3U6  OPT=0 


09/1 c/lZ  lb.2B.07. 


5 


10 


15 


20 


SUBROUTINE  TRAFE 

COMMON/GEOM/RP(6>,X<30) ,R<33) ,C2 , N .NSHAPE, N1 ,N2, XB ( 229) , R3 (225) 
COMMON/GEOl/  RBP  (225)  ,  BET 4 
COMMON/CPU/  CPV<225,20)  ,ja,JO 
COMMON/OIS2/  SUM1,SUM2,SUM3,SUM4 ,SUM5,SUM6 
C  THIS  SUBROUTINE  INTEGRATES  THE  SURFACE  A  RE A , PL ANF ORH  AREA  AND 

C  VOLUME  BY  TRAPEZOIDAL  RULE. 

PI=3. 14159 
IF  (J9.NE. 1)  GO  TO  2 

SUM1=PI*RB(2)*SQRT <RB(2)*»2+XB(?)»*2) 

SUM2=PI/3.*RB(2>*»2+X8<2) 

SUM3=.5*R8(2) *XB (2) 

SUM4=SUH3*2./3.*X3(2) 

GO  TO  99 

2  OC  1  I=JA,JB 

DX=XB ( I  +  1 ) -XP  ( I ) 

SUH1=$U41+PI*0X* (RB(I)*SQRT(l.tRBP(I)**2)tR9(I+l)*SQRT(l.+RPP(I+l 

1 1  **2)  > 

SUM2  =  SUH2*PI/2.*CX» (RB (I) * +2 +RB < I + 1 ) • + 2 ) 

SUM3=SUM3+OX* (RB  (1  +  1)  +RB  (I) ) 

S  UM4=SUM4  +  DX* (XB(I  +  1) * RB ( I +1 >  +XB ( I > *RB ( I > ) 

1  CONTINUE 
99  RETURN 
END 
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SUBROUTINE  NAVE  TRACE  COC  6600  FT N  V3.0-P3U8  OPT=0  09/12/72  16.20.07. 

SUBROUTINE  WAVE 

COMMON/ WAVE/C A  RL  »CNBL  ,CMBL,CAW,CNW,CMK 

COM"ON/GEOM/RP(6> ,X(30),R(30> , C2 , N , NSH APE, N1 , N 2 , X3 ( 225 ) , RB ( 225 ) 

C OMMCN/ GEOl /  RBP(225)  .BETA 

5  COMMON/GE02/NN1,NN2,NN3,NN4,NFL,NBLUNT,NN,NNI,IPRINT,NN1A 

C0MM0N/3IS2/SUM1 ,SUM2 , SUM  3 , SUM4 , SUM  5 , SUM6 
CTMMON/CPV/  CPV(225,20) ,JA,JB 
C0MM0N/GEO4/K,F,RR,RREF 
COMMON/GE03/VOVS.AL , X M , Y M , XI  NT , NN I  A 

10  DIMENSION  F0 (6) ,Fl (6) ,F2 (6) , XN<6> ,RN(6) ,CPN(6,7) ,PH (7) ,RNP (6) , AM (9 

1 )  ,C°0  (O) 

CA2  =  0 . 

CA3=0 . 

CA4=0. 

15  CN2=0 . 

CN3=0 . 

CN4  =  0 • 

CM2=0. 

CM3  =  0 . 

20  C  M4  =  0 ■ 

SUM1=  0 . 

SUM2=  0  • 

SUM3=  0 • 

AREF=3.14159»RREF*»2 

25  IF(NBLUNT.EQ.l)  GO  TO  1 

K-NNI 
K1=NN1 
K2=NN2 
<3=NN3 

30  K4=NN4 

GO  TO  2 
1  K  =  1 

K1=NN1 

K2=NN2 

35  K3=NN3 

K4=NN4 
CABL  =  0 . 

CN8L=0  . 

CM9L=0. 

40  2  JA=K 

J9=K1 
CALL  SIMP 
CA1=  2  .  *SUM1 / ARE F 
CN1=-2.»SUM2/AREF 

45  CM1=  2.*SUM3/(AREF*2.*RRFF) 

SUM  1=  0 • 

SUM2  =  0 • 

SUM3=0  • 

IF1NN1.EC.NN  )  GC  TO  99 
50  J A=K1 ♦ 1 

J  3-K2 
CALL  SIMP 
CA2=  2  •  *  SUM 1 / A  REF 
CN2=-2.'SUME/AREF 

55  CM2=  2.»SUM3/(APEF*2.*KREF) 
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SU9P0UTINT  hAVE  T  pact 


SUM1=  0 • 

SUM2=0 . 

SUM3=0 . 

I F ( NN2 .  EO . NN  )  GC  TO  99 
of  Jfl=K2*l 

J9  =  K3 
can  SIMP 
CA3  =  2 . *  SUM1 /  AREF 
CN3=-2.*SUM2/AREF 

65  CM3=  2.»SUM3/(iOEF*2.»RRC-F) 

SUMisO. 

S UM2=  0  • 

SUM3=  0  . 

IF (NN3.EQ.NN  )  GC  TO  99 
7a  J A=K3 ♦ 1 

J3=K<* 

CALL  SIMP 

cau=  2.*SUMl/aREF 

CN*.  =  -2.*9UM2/APEF 

75  CMi.:  2.*SUM3/<APEF*2.«RREF> 

99  CAW=  CAOLfCai 4CA2FCA3+C49 

CNH=  CN9L*CN1*CN2*-CN3*CN9 
CMW=CMBL»CM1*CM2+CM3+CM4 
PETURN 

80  ENO 


COC  66J 0  FT N  V3.0-P3Q8  OPT=Q  09/12/72  lb.2Q.0T. 


PAGE  2 


E.  SAMPLE  INPUT  DATA  SET 


The  sample  case  below  is  the  input  data  cards  for  the  improved  • 

5"/54  projectile38.  The  Formats  and  parameter  locations  and  definitions 
have  been  discussed  previously. 


F .  SAMPLE  OUTPUT 

The  resulting  output  for  the  above  input  data  case  is  shown  below. 
The  reference  conditions  are  printed  first  followed  by  the  input  body 
coordinates.  The  last  output  quantity  is  a  table  listing  the  indivi¬ 
dual  axial  force  contribution  and  a  table  listing  the  force  coeffici¬ 
ents  . 
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CASE  NO.  1 


ANGLE  OF  ATTACK  s  .500E&S  REFERENCE  01  AN£  TER  s  .417FT 


REFERENCE  CONDITIONS 

SPEED  OF  SOUND  *  1116.890  FT/SEC 
denSITY  *  .0023769  SLUGS/F  T3 

ABSOLUTE  VISCOSITY  «  .000030374528  LB-SEC/FT2 


BODY  COORDINATES 


X 

R 

0.0000 

.0828 

.  J95  0 

.1711 

.790  0 

-.2498 

1.1849 

.3186 

1. 5800 

•  3784 

1.9749 

.4287 

2. 3700 

.4696 

?. 7499 

.5003 

4. 1997 

.5000 

5.2000 

.3683 

*aCH  NO.  S< I N  FRICTION 


?.80f 

.017*. 

2.4P0 

.030? 

2. COn 

.033? 

1.600 

.0360 

1.200 

.0409 

1.050 

•  0426 

1. 000 

.0432 

.95^ 

.0438 

.900 

.0444 

.700 

.043? 

•  500 

.0807 

AXIAL  FORCE  CONTRIBUTIONS 


BASE  PRESSURE 

PRESSURE 

.  0585 

.1063 

.  0656 

.1169 

.  0842 

.1286 

.  1023 

.1415 

.1227 

.1569 

.  1232 

.1527 

.1176 

.0921 

.1037 

•  0304 

.  0666 

.0070 

.  0731 

0.0000 

.  069* 

0.0000 

PROTRUSIONS 

TOJAc 

0.0000 

.  19t«r 

o.caoo 

.2156 

0.0000 

•  246  J 

o.ooco 

.2507 

0. 0080 

.  320  5 

0.0060 

.3195 

0.  000  0 

.  ioi  1 

0.0000 

.1/7  i 

-.0000 

•  13- i 

0.0000 

.120. 

u.OCliO 

.  120- 

N4CM 


N3. 


CO 


CN 


FOPCE  COEFFICIENTS 

CL  CN  CNAL 


chal.  xcp/o 


2.  flOP 

?,uno 

?.cor 

1.60P 

l.POC 
1.050 
1.00C 
.PRO 
.  9C0 
.  YfJO 
.50C 


•  1924 
.2158 
.2462 

•  289  ? 

•  320  5 

•  3 156 
.2530 
.  1750 
.1341 
.1205 
.1205 


.3246 
.0227 
.  0205 
.0121 
•  C131 
.0126 
.0120 
.  0116 
.0115 
.0119 
.0120 


.  0230 
.  0208 
.018*. 
•  Cl  47 

.  0103 

.  00  99 
.  0095 
.  ClOl 
.010^ 
.0109 
.  UllO 


049 

2.823 

-5.o32 

1.9951 

040 

2.605 

-4 .022 

1 .7743 

030 

2.350 

-3. -.58 

1.4717 

017 

1.961 

-1.953 

.9957 

003 

1.496 

-.389 

.2604 

0  04 

1.440 

.494 

-.3409 

009 

1.372 

.978 

-.7120 

Oil 

1.332 

1.301 

-.9772 

010 

1.357 

1.163 

-.8619 

0C2 

1.36b 

-.222 

.1627 

005 

1.376 

-.626 

.3922 
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